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Abstract 


'  Laser  saturation  of  an  atomic  resonance  transition  represents  an 
important  new  mode  of  coupling  laser  energy  into  a  gas  or  plasma.  The 
basic  mechanism  in  either  case  is  superelastic  collisional  heating  of  free 
electrons.  For  a  gas  various  seed  ionization  processes  precede  this 
interaction.  During  the  past  year  we  have  developed  a  computational  code 
for  mapping  the  three  dimensional  nature  of  this  interaction.  This  is 

required  because  this  strong  interaction  invariably  distorts  and  attenuates 
the  laser  pulse  as  it  propagates  through  the  medium  being  excited.  We  have 
also  developed  a  new  experimental  facility  for  studying  this  interaction 
and  have  recently  completed  our  first  spectroscopic  measurements  of  the 
electron  temperature  produced  in  a  sodium  plasma  created  through  laser 
resonance  saturation.  This  temperature  appears  to  be  somewhat  lower  than 
predicted  by  our  computer  simulation  and  we  are  currently  attempting  to 
reconcile  this  difference.  Also  within  the  past  year  we  have  discovered 
that  attenuation  of  the  laser  pulse  is  a  maximum  when  the  laser  is  detuned 

by  about  0.5  nm  from  either  of  the  resonance  lines.  This  result  is  also 

predicted  by  our  computer  code  and  appears  to]  represent  further 
confirmation  of  the  importance  of  laser  energy  extraction  through 

superelastic  electron  heating.  Recently,  we  have  initiated  a  study  of  a 
new  approach  to  the  development  of  a  short  wavelength  laser  based  on  three 
photon  excitation  of  an  ion  with  a  low  lying  metastable  state. 
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RESEARCH  OBJECTIVES 

Laser  saturation  of  an  atomic  resonance  transition  represents  an 
Important  form  of  laser-plasma  interaction.  Even  In  the  case  of  a  neutral 
atomic  vapor  laser  resonance  saturation  gives  rise  to  many  processes  that 
quickly  lead  to  ionization.  Once  an  appreciable  degree  of  ionization  has 
been  achieved  the  primary  mechanism  for  converting  laser  energy  into  plasma 
energy  is  superelastic  electron  collision  quenching  of  the  laser  pumped 
resonance  state  population.  This  is  an  extremely  efficient  process  which, 
in  effect,  rapidly  cycles  the  bound  electron  between  the  ground  and 
resonance  state,  transforming  one  laser  photon  into  an  equivalent  amount  of 
free  electron  kinetic  energy  during  each  cycle,  see  Fig.  1. 

Laser  resonance  saturation  can  be  viewed  as  complementary  to  inverse 
bremsstrahl ung ,  for  although  it  can  lead  to  energy  deposition  that  is  many 
orders  of  magnitude  greater  than  attainable  through  inverse  bremsstrahl ung 
for  modest  values  of  laser  irradiance  —  the  free  electron  temperature 
achieved  is  limited  by  the  characteristics  of  the  atom  or  ion  being 
irradiated.  However,  laser  resonance  saturation  can  lead  to  extremely 
rapid  electron  heating,  in  effect,  creating  a  well  defined  electron 
temperature  jump  and  a  commensurate  sharp  burst  in  the  emission  of  short 
wavelength  radiation  from  a  suitable  plasma. 

The  objective  of  our  research  progam  is  to  study  this  new  form  of 
laser-plasma  interaction  and  explore  some  of  its  potential  applications. 
In  any  practical  application  the  laser  has  to  propagate  through  a  finite 
depth  of  the  medium  being  irradiated  and  as  a  result  of  the  strong  nature 
of  the  resonance  interaction  the  laser  pulse  will  undergo  an  appreciable 
degree  of  attenuation  and  temporal  and  spatial  distortion.  Consequently, 
any  theoretical  modeling  of  this  interaction  should  take  account  of  this 
fact. 

To  date  we  are  the  first  to  develop  a  computational  analysis  that 
attempts  to  model  the  three  dimensional  nature  of  this  interaction  by 
taking  account  of  the  attenuation  and  distortion  suffered  by  the  laser 
pulse  as  it  burns  its  way  through  a  realistic  atomic  vapor  distribution. 
In  order  to  assess  the  reliability  of  our  theoretical  work  we  are  comparing 
our  computer  simulations  with  the  results  attained  from  experiments 
involving  a  new  facililty. 
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Multiphoton  saturation  of  atomic  transitions  within  a  plasma 
represents  the  logical  extension  of  this  work  and  we  have  initiated  a  new 
program  to  consider  one  possible  application  of  this  interaction.  This 
work  Involves  examining  the  feasibility  of  a  new  approach  to  the 
development  of  a  short  wavelength  laser  based  on  three  photon  saturation  of 
an  ion  containing  a  low  lying  metastable  level. 

STATUS  REPORT 

Lasers  have  been  used  to  create  and  heat  plasma  almost  since  their 
inception.  In  the  past  inverse  bremsstrahl ung  has  represented  the  primary 
means  of  heating  the  free  electrons  of  a  plasma.  Unfortunately,  the  small 
cross  section  involved  in  this  interaction  has  limited  laser  heating  to 
high  density  plasmas.  By  contrast  laser  resonance  saturation  involves  a 
very  large  cross  section  so  that  efficient  coupling  of  laser  energy  can  be 
achieved  at  much  lower  densities.  In  this  section  we  shall  report  on  our 
work  to  study  laser  resonance  saturation  and  explore  its  uses. 

Laser  Ionization  Based  on  Resonance  Saturation  (LIBORS) 

An  intense  laser  field  of  appropriate  wavelength  can  lock  the 
populations  of  the  ground  and  resonance  states  in  the  ratio  of  their 
degeneracies  in  an  extremely  short  time.  This  laser  saturation  vrtien 
applied  to  an  atomic  resonance  transition  has  been  shown  to  produce 

near  full  ionization  and  there  is  now  substantial  support  for  our  idea  that 
the  primary  mechanism  for  coupling  the  laser  radiation  to  the  free 
electrons  is  superelasic  collision  quenching  of  the  laser  pumped  resonance 
state.  As  first  pointed  out  by  the  author^'*^  the  resonance 

transition,  in  effect,  behaves  as  an  extremely  efficient  engine  that 
converts  laser  energy  into  kinetic  energy  of  the  free  electrons,  see  Fig. 

1. 

We  have  shown  theoretically^"^  that  this  new  form  of  laser-plasma 
interaction  should  lead  to  a  jump  in  the  free  electron  temperature  that  is 
independent  of  laser  irradiance  (provided  it  Is  adequate  to  maintain 
resonance  saturation)  and  limited  in  magnitude  by  the  characteristics  of 
the  atom  or  ion  being  excited.  The  risetime  for  the  temperature  is 
principally  determined  by  the  density  of  the  plasma.  An  attractive  feature 
of  superelastic  heating  is  the  efficient  production  of  a  concomitant  flash 
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of  intense  line  radiation.  This  high  efficiency  is  inherent  in  the 
mechanism  that  limits  the  rise  In  the  free  electron  temperature  —  that 
is  to  say  a  dynamic  balance  between  electron  superelastic  and  inelastic 
collisions.  The  former  heat  the  free  electrons  while  the  latter  excite 
higher  energy  levels  at  the  expense  of  energy  from  the  free  electrons.  In 
effect,  the  laser  energy  is  being  converted  into  intense  line  radiation 
with  high  efficiency. 

The  wavelength  of  this  radiation  can  be  very  much  shorter  than  the 

wavelength  the  saturating  laser  if  the  selected  ion  has  an  appropriate 

energy  level  structure.  Members  of  the  Li-isoelectronic  series  (Bell, 

Bill,  CIV,  NV,  0V1 ,  ...)  are  near  ideal  as  they  possess  a  large  energy  gap 

between  their  first  and  second  resonance  levels.*  Table  1  displays  the 

wavelengths  of  the  two  lowest  resonance  transitions  for  several  members  of 

this  series  indicating  clearly  that  the  wavelength  of  the  intense  line 

emission  would  be  very  much  shorter  than  that  of  the  saturating  laser. 

This  photon  energy  gain  factor  (x,,/X3])  can  be  seen  to  increase 

/pi  i 

progressively  with  the  members  of  the  series.  We  have  also  shown'  '  that 
the  magnitude  of  the  temperature  jump  achieved  through  resonance  saturation 
increases  in  a  similar  manner. 


Table  1 

Resonance  Lines  of  the  Li-I soelectronic  Series 
and  Photon  Energy  Gain  Factor  (\2iA3i) 


Bell 

Bill 

CIV 

NV 

0VI 

X31 

103.6  nm 

51.8  nm 

31. 2  nm 

20. 9  nm 

15. 0  nm 

*21 

313  nm 

207  nm 

155  nm 

124  nm 

103  nm 

X  2 1  /  X  3 1 

3.0 

4.0 

4.97 

5.93 

6.87 

It  is  worth 

noting  at 

this  point 

that  the 

sel ection 

of  wavelengths 

attainable  through  superelastic  heating  is  not  restricted  to  those 
*  We  are  neglecting  the  fine  structure  of  the  levels  in  this  statement. 


associated  with  the  members  of  this  series  for  they  could  be  used  as  a 

source  of  temperature  jump  within  a  plasma  and  another  species  serve  to 

provide  the  specific  line  radiation  required. 

We  believe  that  this  new  kind  of  efficient  and  controlled  "intense 

short  wavelength  line  emission"  (ISLE )-source  could  be  used  in  many  ways 

towards  the  development  of  short  wavelength  (EUV  to  X-ray)  lasers.  For 

(22) 

example,  in  the  approach  of  Caro  et  al '  '  an  ISLE-source  that  is  tailored 

to  provide  primarily  radiation  at  18.7  nm  (as  opposed  to  the  broadband 

source  they  used)  would  avoid  excitation  of  the  Li  states  that  subsequently 

autoionize  to  the  Li+  ground  state  and  destroy  the  chance  of  creating  a 

population  inversion.  A  suitably  chosen  ISLE-source  would  also  be  ideal 

( 

for  the  selective  autoionization  proposal  of  Bokor  et  al.  ' 


Three  Dimensional  Modeling  of  L I BOR S 

In  order  to  better  understand  laser  ionization  based  on  resonance 
saturation  we  have  undertaken  the  first  theoretical  treatment  that  attempts 
to  model  the  three  dimensional  nature  of  the  interaction  and  includes  such 
real  effects  as:  laser  pulse  attenuation,  non-uniform  density  distribution 


along  the  path  of  the  laser  pulse  and  spectral  detuning  of  the  laser. 

In  our  original  formulation  of  the  L1B0RS  computer  code  we  assumed  a 
step-like  laser  pul  se.  This  approximation  is  reasonable  for  the 
incident  laser  pulse  because  the  time  to  saturate  the  resonance  transition 
is  very  short  typically  much  less  than  1  nsec.  However,  in  an 


actual  experiment  the  laser  pulse  propagates  through  a  non-uniform  atom 
density  distribution  and  suffers  appreciable  attenuation.  This  can  lead  to 


a  consioerable  temporal  and  spatial  distortion  of  the  laser  pulse. 
Consequently,  the  atoms  at  different  locations  along  the  path  of  the  laser 
pulse  are  subject  to  progressively  less  energy  and  to  a  shorter  period  of 


excitation. 


5-Level  LIBORS  Computer  Code 

In  an  attempt  to  keep  the  analysis  economically  viable  we  have 
replaced  our  20-level  model  0f  the  sodium  atom  with  the  simplified 
5-level  model  shown  in  figure  2.  Although  5-levels  are  indicated  in  figure 
2,  we  assume  that  the  collisional  coupling  between  levels  3,  4  and  5  is 
sufficiently  rapid  that  their  populations  are  more  or  less  in  a  Boltzmann 
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equilibrium.  This  assumption  permits  us  to  treat  this  trio  of  levels  as  a 
manifold  and  thereby  reduce  the  number  of  atomic  level  rate  equations  to 
three.  In  effect  we  consider _only  the  collisional-radiative  processes 
"into"  and  "out  of"  this  (3-4-5)  manifold. 

The  remaining  15  high  lying  energy  levels  (not  shown  in  figure  2)  are 
partially  taken  into  account  by  assuming  that  electron  excitation  to  them 
from  any  of  the  five  lower  levels  leads,  in  effect,  to  immediate  ionization 
of  the  atom.  We  therefore  introduce  the  5-level  "effective"  electron 
collisional  ionization  rate  coefficient  for  each  of  the  low  lying  levels  (m 
=  1  to  5)  defined  by  the  relation 


K* 

me 


21 

l  K 

n=6  mn 


(1) 


where  Kmn  represents  the  electron  collisional  rate  coefficient  between 
levels  m  and  n,  as  used  in  our  LIBORS  computer  code.  Note  that 

is  the  m-level  electron  ionization  rate  coefficient  according  to  the 
20-level  model  of  the  sodium  atom.  The  corresponding  "effective"  three 
body  recombination  rate  coefficient  to  each  of  the  low  lying  levels  (m  =  1 
to  5)  is  similarly  constituted  to  take  account  of  recombination  to  the 
higher  levels,  followed  by  a  rapid  cascade  down  to  the  lower  levels,  we 
thus  write 


20 

K*  =  T  K  p(n,m)  (2) 

cm  Lc  cnK'  ’  '  v  ' 

n=o 

where  Kcn  represents  the  three-body  recombination  rate  coefficient  to  level 
n  used  in  our  20-level  LIBORS  code  and  p(n,m)  is  introduced  as  the 
effective  probability  that  a  captured  free  electron  arriving  in  level  n 
subsequently  decays  to  level  m. 

We  have  found  that  close  agreement  can  be  obtained  between  the  5  and 
20-level  LIBORS  codes  if  these  probability  factors  are  generally*  given  by 
the  expression 

*  In  order  to  obtain  the  close  agreement  (indicated  later)  between  the  5 
and  20-level  LIBORS  codes  it  was  found  necessary  to  set  p(s,  4)  =  1,  where 
s  =  17,  12  and  8. 
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(3) 


■ 

E 


g 


* 


p(n,m) 


IT  (If 
e '  nm 


N  K  ♦  A 
e  nm  nm 

TT 


„„  tt  y~rj  +~A 
nq  nc '  nm  nq 


where  Is  the  free  el ectron- densi ty ,  Anm  the  Einstein  spontaneous 
transition  probability  between  levels  n  and  m,  and  q  is  the  level  within 
the  manifold  that  is  optically  connected  to  level  n. 

The  5-level  model  can  thus  be  described  by  the  following  set  of 
equations:  first  we  have  the  three  atomic  level  population  rate 

equations , 

dN. 

=  (Nj  -  gffj)R21  *  N?C2j  *  NjCy  *  Ne(N3K31  *  N4K41  )  *  Ne2Pcl 

4NAWNen  Klm  *  Kiel]  <4> 

m=6 


dN, 


ar  *  <9N1  -  N2>R21  *  Wl2  4  N5K52>  *  N3C32  4  N4C42  4  Ne2f,c2 


*  Ne3Kc*2  -  W,  l  *2™  *  K2*c>  4  C21 )  ‘  7  N2^<1F  *  »A>  ‘  N2<4I>F2 
m= 6 


(5) 


and 


dN3 


■ari  *  N  [  y  (N.K,  *  N,K,  *  N  V  *  N  0  )]  -  N, (C,,  +  N  K,,  ♦  P,  ) 

dt  e  m=3  1  2  2m  e  cm  e  cm,J  3 v  32  e  31  3c7 


-  VC42  4  NeK41  4  P4c>  -  VC51  4  "eK52  4  P6c>  <6> 


where  we  have  introduced  the  (3-4-5)  level  manifold  population  density 


N*  =  N,  +  N,  +  N, 
3  3  4  5 


(7) 
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4 


r 


As  Indicated  we  assume  that  the  close  spacing  of  these  three  levels  ensures 
collisional  equilibrium  and  consequently  that  the  internal  distribution  of 
this  manifold  is  described  by  a  Boltzmann  distribution.  In  which  case 


where 


Nk  /Uk(Te) 

-WkTe 


Uk(Te)  *  moJ.S  U' 


(8) 

(9) 


with  =  gm/gk  and  k  =  3,  4  or  5.  In  order  to  streamline  the  equations 

as  far  as  possible  we  also  introduced  the  notation, 


and 


nm 


■  NeK™ 


+  A. 


nm 


(10) 


nc 


=  o  F  +  N  K* 
nc  e  nc 


(11) 


where  onc  is  the  single  photon  ionization  cross  section  for  level  n  and  F 
is  the  laser  photon  flux  density,  i.e., 

F(t)  -  -ji  /  lV,t)dv  (12) 

where  I  (v,t)  is  the  laser  spectral  irradiance  at  time  t  and  hv  the  laser 
photon  energy.  In  the  above  equations  we  have  also  introduced  the 
resonance  to  ground  level  degeneracy  ratio  g  and  the  stimulated  emission 
rate  per  resonance  state  atom, 

821  , 

^21  (t)  =  /  1  { v ,t ) L g j  (v)dv  (13) 

where  is  the  Milne  stimulated  emission  coefficient  and  (v)  is  the 
line  profile  function  for  the  resonance  transition. 

The  appropriate  electron  energy  equation  can  be  expressed  in  the 
following  form: 
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O 

where  eg  represents  the  free  electron  mean  thermal  energy  (taken  as  kTg), 

Pcn  is  the  radiative  recombination  rate  coefficient  appropriate  for  level 

n»  H  .  represents  the  rate  of  elastic  energy  transfer  to  ions  through 
ei  (?&\ 

Coulomb  scattering  collisions'  '  and  S  corresponds  to  the  source  of 
electron  energy  arising  from  the  various  seed  ionization  processes,  i.e.. 


S  =  (2E21  *  Ec2 ^ ^2°2c)p 2  +  7  U22yai^  +  7  N22yoAEA  (15) 

Here  o£)  represents  the  two  photon  resonance  state  ionization  coefficient, 
oL  the  laser  induced  Penning  ionization  coefficient,  the  associative 
ionization  cross  section  and  E^  the  corresponding  energy  of  creation  for 
the  electrons  liberated  by  this  latter  process,  v  represents  the  mean  atom 
vel ocity . 

The  energy  equation  is  subsequently  divided  into  an  ionization 
equation  and  an  electron  temperature  equation  as  discussed  in  our  earlier 
work.^17-20^  It  should  be  noted  that  another  fine  adjustment  in  the  form 
of  a  small  reduction  (of  about  0.055  eV)  in  the  three-body  recombination 
energy  was  found  to  bring  the  results  of  the  5-level  LIBORS  code  into 
excellent  agreement  with  the  20-level  code  over  a  wide  range  of  densities. 
This  modified  energy  is  designated  E*  in  equation  (14). 

As  we  shall  see  shortly  the  resonance  state  population  density  N2  is 
critical  in  establishing  the  extent  of  laser  attenuation  and  its  value  is 
in  turn  determined  by  the  free  electron  temperature  Tg.  Once  appreciable 
ionization  has  been  attained  Tg  tends  to  rise  above  the  low  value  produced 
as  a  result  of  burnout.  This  is  particularly  true  at  the  higher  densities 
(Nq  >  1016  cm-2).  The  electron-ion  exchange  term  is  the  primary  limit 
to  the  growth  in  T  .  Since  is  proportional  to  (Tg-T^ )  it  is  necessary 
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to  also  Introduce  the  appropriate  ion  energy  equation 


W-  \  NeHei  ‘  NAH1A  (16) 

where  Is  the  mean  thermal  energy  of  the  Ions  (set  equal  to  kT^  ) .  The 
second  term  In  equation  (16)  takes  account  of  elastic  cooling  collisions 
between  the  sodium  ions  and  the  background  of  argon  buffer  gas  atoms.  NA 
is  the  background  argon  atom  density  and  H^A  is  the  rate  of  elastic  energy 
exchange  between  the  Ions  and  the  background  argon  atoms,  i.e.. 


HiA  =  «rv>< 


4m.. 


iA  m, 


-ta> 


A  i 


(17) 


where  <ov>iA  is  the  ion-argon  collision  frequency  and  TA  the  argon  atom's 
translational  temperature. 

This  background  density  of  argon  atoms  varies  across  the  sodium  vapor 

region  such  that  the  sum  of  the  partial  pressures  equals  the  pressure  of 

the  argon  in  the  water  cooled  periphery  region.  The  initial  vapor 

temperature  at  each  location  is  deduced  from  the  optically  measured  sodium 

f  2  5 1 

density  profile  using  the  vapor  pressure  curves  provided  by  Nesmeyanov.  ‘ 
By  a  similar  argument  to  that  given  above  it  follows  that  since  HiA  is 
proportional  to  the  difference  in  temperature  between  the  sodium  ions  and 
the  aroon  atoms,  the  latter  has  also  to  be  computed.  This  is  done  by  means 

of  an  argon  energy  equation.  It  should  also  be  noted  that  an 

electron-argon  elastic  energy  loss  term  was  also  introduced  into  both  the 

free  electron  energy  equation  and  the  argon  energy  equation  but  its 

relative  importance  is  marginal. 

Combining  equations  (12)  and  (13)  enables  us  to  write 


R2i(z»t)  F(z»t;)  /  LA(xA-x)L2i(x)dx  (18) 

where  X  is  the  wavelength  and  x21  the  radiative  lifetime  of  the  resonance 
transition.  We  have  also  introduced  the  spectral  distribution  function  for 
the  laser  radiation 
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(19) 


■  1 1  .  .  IA(v) 

L  <VV>  =-TuT 


and  the  dummy  frequency  variable  _ 


x  =  v  -  v. 


where  vq  is  the  resonance  line  centre  frequency.  For  our  experiment  we  can 
assune  a  Lorentzian  line  profile  function  of  the  form 


L2l(v  -  v  J  «  - 


(v-vj2  +  Y, 


where  y  is  the  appropriate  HMHW  bandwidth  (the  Stark  width  in  essence). 
The  laser  spectral  distribution  function  can  well  be  approximated  by  a 
Gaussian 


and  consequently  the  convolution  of  the  two  profiles  in  (18)  can  be 
expressed  in  the  form  of  a  Voigt  function, 


V(o,u) 


e'^dv 


a2  +  (u-y) 2 


where 


a  =  ya/0(,  y  =  (xA-x) u  =  Xf/P,  and  V(a,  u)du  =  V(a,  v)dv 


It  is  important  to  notice  that  this  Voigt  function  depends  upon  a  -  the 
ratio  of  the  atomic  bandwidth  to  the  laser  bandwidth  and  to  the  frequency 
separation  of  the  atomic  and  laser  centre  frequencies  normalized  by  the 
laser  bandwidth,  i .e. , 
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(24) 


In  the  early  phase  of  our  -work  we  were  interested  in  comparing  our 
theoretical  predictions  with  experiments  based  on  the  use  of  a  Phase-R 
flashlamp  pumped  dye  laser,  the  output  of  which  had  a  duration  of  about  1 
usee.  A  good  fit  to  the  temporal  laser  pulse  from  this  system  was  obtained 
using  the  empirical  expression 


F(o,t)  =  Fmax(t/400)l'2exp{-  j  (t/400)  (25) 


where  Fmflx  is  the  peak  laser  photon  flux  density  and  t  is  in  nanoseconds. 
The  exponential  bandwidth  of  this  laser  pulse  was  measured  to  be  about 
1011  Hz  (corresponding  to  a  full  exponential  linewidth  of  around  0.25  nm). 

We  recognized  from  the  outset  that  attenuation  of  the  laser  pulse 
could  lead  to  relatively  sudden  changes  in  the  laser  irradiance. 
Consequently  in  our  comparison  of  the  5  and  20-level  computer  codes  we  used 
the  truncated  laser  pulse  indicated  as  F  in  figure  3(a).  With  a  laser 
pulse  of  this  shape  the  comparison  was  tested  fairly  severely. 

The  results  for  a  sodium  density  (Nq)  of  1016  cm* 3  are  presented  in 
figure  3(a,  b  and  c).  It  is  evident  from  a  comparison  of  figures  3(a)  and 
3(b)  that  N2/N 2  falls  from  the  saturated  value  of  3.0  at  about  the  time  of 
ionization  burnout  even  though  F  is  still  increasing.  This  temporal 
minimum  is  attributed  to  the  sudden  increase  in  the  superelastic  quenching 
rate.  It  is  apparent  from  figure  3  that  the  5-level  LIBORS  code  mirrors 
almost  exactly  the  predictions  of  the  20-level  code  in  regard  to  the 
temporal  behaviour  of  the  various  parameters  (i.e.,  Ng,  T  ,  T^,  and  Nj). 
We  feel  that  the  small  (*  20%)  difference  in  the  lead  time  between  the  two 
sets  of  curves  is  quite  acceptable.  In  order  to  check  that  this  close 
agreement  was  not  fortuitous  to  this  particular  density  we  ran  a  comparison 
at  two  additional  densities  (NQ  =  5*1015  and  4*1016  cm*3).  The  higher 
density  comparison  was  also  run  for  a  much  shorter  laser  pulse, 
corresponding  to  our  new  laser  system,  and  is  presented  as  figure  4.  It  is 
evident  that  the  same  close  agreement  is  achieved.  Similar  close  agreement 
was  obtained  for  the  lower  density  comparison. 
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The  close  agreement  between  the  two  codes,  even  vrfien  the  laser 
irradiance  was  suddenly  truncated,  provided  us  with  sufficient  confidence 
to  use  the  5-level  LIBORS  code  In  evaluating  the  time  history  of  Ionization 
and  the  attenuation  of  the  laser  pulse  as  It  propagates  through  the 
extended  region  of  sodiun  vapor. 

Laser  Pulse  Attenuation  Within  an  Extended  Atomic  Vapor  Reoion 


If  we  treat  the  problem  as  one  dimensional  the  radiative  transfer 
equation  governing  the  propagation  of  the  laser  pulse  can  be  written  In  the 
form 

*  £  ■  <N2  ’  9N1>R21  -  V2c,f2  -  7  "2V  *  l  W  <26> 

m= o 

The  first  term  on  the  right  hand  side  takes  account  of  the  loss  of  laser 
photons  due  to  absorption  in  the  resonance  transition,  the  second  term 
allows  for  two  photon  ionization  of  the  resonance  level,  the  third  term 
involves  laser  induced  Penning  ionization  of  the  resonance  level  and  the 
last  term  refers  to  single  photon  ionization  of  the  levels  within  the 
mani fold. 

For  the  time  scale  of  interest  in  our  experiments  the  velocity  of 
light  c  can  be  treated  as  infinite.  Under  these  circumstances  (26) 
is  simplified  to  the  steady  state  form  given  by 

£  ■  ("2  '  SN1>R21  ‘  N2°2c>fZ  ’  7  VV  -  \  NmVF  <”> 

m=  j 


where  R?^  was  given  by  equation  (18). 

It  is  quite  clear  that  addressing  the  problem  of  laser  resonance 
saturation  for  an  extended  vapor  region  involves  the  solution  of  equations: 
4,  5,  6,  14,  16  and  27  in  both  space  and  time.  We  decided  to  tackle  this 
problem  by  dividing  the  vapor  column  along  the  path  of  the  laser  pulse  into 
a  series  of  slices,  see  figure  5,  each  of  which  was  assumed  to  have  a 
uniform  atom  density  equal  to  the  average  experimentally  measured  sodium 
atom  density  for  that  segment  of  vapor.  The  thickness  of  each  slice  was 


initially  set  at  1  on  although  later  computer  runs  used  a  much  finer  slice 
di strlbution. 

Within  each  slice  our  5-level  L1B0RS  computer  code  was  used  to 
evaluate  the  temporal  variation -of  the  various  parameters  (such  as  N  ,  T  , 
Tj »  NJ,  N2  and  Nj )  subject  to  the  laser  pulse  exiting  the  previous  vapor 
slice.  Although  attenuation  of  the  laser  pulse  within  each  slice  is 
neglected  [i.e.,  we  treat  F(z,t)  as  F(t)],  the  energy  extracted  within  this 
vapor  slice  is  computed  from  equation  (27)  and  this  loss  Is  taken  into 
account  when  considering  the  appropriate  laser  pulse  to  use  for  the 
subsequent  vapor  si  ice. 

Representative  time  histories  of  the  free  electron  density,  the 
electron  temperature  and  the  laser  pulse  photon  flux  density  at  the  three 
locations  indicated  on  the  sodium  atom  density  distribution  (measured  in  an 
actual  experiment  and  displayed  in  the  lower  left-hand  corner)  are 
presented  in  figure  6.  The  incident  laser  pulse  used  in  this  simulation 
was  chosen  to  closely  match  that  obtained  from  our  new  Nd-YAG  laser  pumped 
dye  laser.  An  excellent  fit  to  this  short  laser  pulse  was  provided  by  the 
empirical  expression 

F(o,  t)  -  F  x  1.79  x  10-2  t2*06  exp(-7.  9*10" 3t2)  (28) 

where  t  is  in  nanoseconds,  see  figure  7.  For  the  computer  simulation 
presented  in  figure  6,  Fp  =  8. 33xl025  photons  cm" 2  sec"1,  corresponding  to 
an  energy  of  60  mJ  in  0.157  cm2  beam  cross  section,  and  the  laser  was 
assumed  to  be  detuned  by  0.14  nm  from  the  589.0  nm  resonance  line. 

It  is  quite  apparent  that  at  this  energy  full  ionization  is  predicted 
across  the  entire  vapor  column  and  that  the  general  shape  of  the 
temperature  profile  should  be  quite  similar  at  each  location  —  albeit 
lengthened  or  shortened  in  time  dependintg  upon  the  location.  The  time  to 
achieve  full  ionization  is  also  seen  to  vary  across  the  vapor  column,  being 
a  minlmun  where  the  density  is  a  maximum. 

Unfortunately,  the  smaller  the  spectral  detuning  of  the  laser  from  one 
of  the  sodium  resonance  lines  the  stronger  the  radiative  rates  R12  and  R2i 
and  the  smaller  the  time  and  space  step  size  required  for  the  computer 
calculations.  Experience  has  Indicated  that  the  smallest  economically 
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acceptable  value  of  laser  detuning  was  0.14  nm  for  conditions  commensurate 
with  our  experiments. 

The  spatial  variation  of  the  degree  of  ionization  along  the  vapor 
column,  irradiated  by  the  laser  pulse,  predicted  at  several  instants  from 
the  onset  time  of  the  incident  laser  pulse  (of  energy  27  mJ  and  detuned  by 
0.14  nm  from  the  589.0  nm  sodium  resonance  line)  is  displayed  in  figure  8. 
It  is  evident  that  even  with  this  reduced  energy  near  complete  ionization 
is  achieved  over  most  of  the  vapor  column  by  the  end  of  the  laser  pulse  (30 
ns). 

In  figure  9,  the  spatial  variation  of  the  electron  temperature  and  the 
degree  of  ionization  predicted  to  exist  along  the  irradiated  vapor  column 
10  ns  from  the  start  of  the  incident  laser  pulse  for  three  values  of  laser 
detuning  (0.14,  0.4  and  2.0  nm)  are  presented.  Large  gradients  of  both 
electron  temperature  and  density  are  clearly  predicted  for  laser  detunings 
of  0.14  nm  (or  less).  Indeed,  it  is  expected  that  steep  fronts  of 
ionization  and  electron  temperature  propagate  along  the  irradiated  vapor 
column,  especially  for  zero  detuning  of  the  laser. 

Unfortunately,  at  about  this  time  the  laser  had  to  be  realigned  after 
a  component  failure  and  this  realignment  caused  the  laser  pulse  to  become 
slightly  longer  in  duration  (40  ns  compared  to  30  ns)  and  larger  in  cross 
section.  The  new  empirical  fit  to  the  laser  temporal  profile  was  found  to 
be  given  by 

F (o,  t)  =  FpAtBe'Ct2  (29) 

where  A  =  3.352xl0"3,  B  =  2.642  and  C  =  6.50xl0"3.  The  new  cross  sectional 

area  was  0.236  cm2  and  the  current  relation  between  the  peak  photon  flux 

£ 

density  Fp  and  the  laser  energy  fluence  e  is  of  the  form 

Fp  (photons  cm"2  sec"1)  *  2xl023  e*  (mJ  cm"2)  (30) 

In  order  to  model  the  three  dimensional  nature  of  this  interaction  we 
divided  the  lateral  distribution  of  the  laser  pulse  into  several  concentric 
annular  shells  and  equated  the  laser  energy  fluence  of  each  shell  to  the 
corresponding  mean  value  in  the  radial  distribution. 


The  Interaction  along  the  length  of  each  shell  was  then  computed.  The 
lateral  variation  of  electron  density,  electron  temperature  and  laser 
Irradlance  could  then  be  determined. 

An  example  of  the  lateral-  distribution  of  electron  density  and 
electron  temeprature  predicted  at  the  centre  of  the  sodium  vapor 
distribution  for  two  instants  of  time  Is  displayed  in  figure  10.  In  this 
instance  the  computer  was  simulating  an  experimental  run  in  which  the 
sodium  atom  density  at  the  centre  of  the  heat  sandwich  oven  was  measured  to 
be  1.5*1016  cm' 3  and  the  energy  in  the  laser  pulse  was  estimated  to  be  19.7 
mJ  with  a  burn  spot  average  radius  of  2.74  rrm.  For  the  computer  run  the 
incident  laser  pulse  was  assumed  to  have  a  Gaussian  radial  distribution 
with  an  exponential  radius,  rQ  of  1.225  mm.  The  incident  peak  radial  laser 
photon  flux  density  F*  distribution  is  also  shown  in  figure  10. 

It  is  apparent  that  the  computer  code  predicts  that  near  full 
ionization  is  achieved  out  to  a  radius  of  about  1.5  mn  and  that  the 
electrons  are  hottest  at  the  edge  of  the  laser  beam.  The  predicted  cool 
annular  region  adjacent  to  the  hot  region  is  associated  with  a  rapid 
ionization  burnout  occurring  at  25  ns  at  that  radial  location.  It  should 
be  mentioned  that  our  assumption  of  a  Maxwellian  distribution  for  the  free 
electrons  is  likely  to  break  down  for  N&  <  1015  cm'3  and  so  the  reliability 
of  our  temperature  predictions  very  close  to  the  edge  of  the  laser  beam  is 
questionable. 

As  indicated  in  the  subsequent  section  we  are  on  the  verge  of 
measuring  the  radial  electron  density  and  temperature  profiles  produced  by 
laser  resonance  absorption.  Consequently,  as  an  initial  test  of  the 
accuracy  of  our  three  dimensional  modeling  of  this  laser  ionization 
technique  we  decided  to  compare  the  predicted  temporal  shape  and  energy 
fluence  of  the  transmitted  laser  pulse  with  that  observed  in  an  experiment. 
From  the  many  computer  runs  corresponding  to  various  incident  laser  energy 
pulses  we  have  drawn  a  curve  of  laser  energy  fluence  absorbed  (nvJ  cm-2) 
against  incident  peak  laser  photon  flux  density  Fp  (or  energy  fluence)  for 
the  case  of  a  peak  sodium  density  of  1. 5*10 1 6  cm" 3  and  a  laser  detuned  to 
lie  halfway  between  the  two  resonance  line  wavelengths,  i.e.,  X  =  589.3 
nm.  This  curve  is  presented  as  figure  11,  and  using  it  we  can  calculate 
the  amount  of  laser  energy  transmitted  in  each  of  the  annular  shells.  This 


enables  us  to  essentially  predict  the  amount  of  energy  transmitted  through 
the  above  sodium  atom  distribution  for  any  radial  profile  of  incident  laser 
energy  of  wavelength  589.3  nm. 

To  compute  the  temporal  shape  and  energy  of  any  given  transmitted 
laser  pulse  we  must  reconstitute  this  laser  pulse  from  the  set  of 
transmitted  laser  pulses  that  were  used  to  represent  the  incident  laser 
pulse.  The  results  of  our  first  such  synthesis  and  the  corresponding 
experimental  data  are  presented  in  figure  12.  It  can  be  seen  that  the 
predicted  transmitted  laser  pulse  approximates  quite  well  the 
experimentally  monitored  laser  pulse  propagating  through  the  sodium  vapor. 
This  preliminary  comparison  gives  us  some  degree  of  confidence  in  our  three 
dimensional  modeling  of  LIBORS. 

New  LIBORS  Facility  and  Experiments 

During  the  past  year  our  laboratory  has  undergone  a  major 
reorganization  centred  about  the  installation  of  a  new  laser  facility  and 
an  associated  electronic  data  acquisition  and  reduction  system.  As  stated 
in  our  previous  technical  report  the  Phase-R  flashlamp  pumped  dye  laser 
which  had  formed  the  workhorse  of  our  initial  LIBORS  experimental  program 
suffered  several  major  failures  and  never  met  specifications. 

At  the  heart  of  the  new  facility  is  a  Nd-YAG  HY750 
oscillator/ampl ifier  laser  manufactured  by  J.  K.  Lasers  of  the  U.K.  and 
distributed  on  this  continent  by  Lumonics.  This  system  also  includes 
second  and  third  harmonic  generators  and  its  specifications  are  indicated 
in  Table  2.  The  second  part  of  this  faclity  is  a  Quanta-Ray  (oscillator/ 
amplifier)  PDL-1  dye  laser  that  is  pumped  by  either  the  second  or  third 
harmonic  of  the  Nd-YAG  laser.  We  also  have  provision  for  frequency 
doubling  the  output  of  the  dye  laser  or  frequency  mixing  it  with  the 
fundamental  of  the  Nd-YAG  laser.  An  overview  of  this  new  laser  facility  is 
presented  as  figure  13,  while  a  photograph  of  the  laboratory  is  shown  as 
figure  14.  The  Nd-YAG  laser  can  be  seen  in  the  background.  The  sodium 
heat  sandwich  oven  with  its  flexible  tubing  carrying  the  cooling  water  is 
clearly  visible  at  the  centre  of  the  picture.  The  photodetection  system, 
including  the  vertically  mounted  SPE X  17001 1  monochromator,  a  small  Heath 
monochromator,  various  photodiodes  and  a  camera,  are  also  seen  on  the 
granite  table. 
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Output  energy 

Nd-YAG 

(at  1.06  pm) 

850 

mJ 

(at  532  nm) 

320 

mJ 

(at  355  nm) 

220 

mJ 

FWHM  pulse  duration 

(at  1.06  pm) 

20 

ns 

(at  532  nm) 

17 

ns 

Repetition  rate 

10 

Hz 

In  order  to  optimize  the  use  of  this  versatile  and  powerful  new  laser 
facility  we  have  also  switched  our  mode  of  operation  and  data  acquisition 
during  this  period.  Originally  we  recorded  single  pulses  on  a  fast 
oscilloscope  using  cameras,  then  manually  digitized  these  photographs. 
With  the  purchase  of  an  EGG  (4420)  signal  averager  we  now  acquire  the 
photomultiplier  signals  electronically  and  perform  data  reduction  by  means 
of  an  Apple  II  Plus  microcomputer  that  also  serves  as  a  terminal  to  the 
Institute's  new  Perkin  Elmer  3250  computer.  This  reorganization  of  the 
laboratory  has  limited  our  experimental  program  during  the  first  six 
months,  but  in  the  past  six  months  we  have  completed  several  runs  with  this 
new  facility. 

An  evaluation  of  the  free  electron  temperature  of  the  sodium  plasma 
produced  through  laser  resonance  saturation  constitutes  one  of  the  most 
important  measurements  to  be  undertaken.  Initially  the  emission  from  a 
series  of  spectral  lines  was  monitored  and  the  spatially  averaged  electron 
temperature  deduced.  This  measurement  was  undertaken  by  imaging  a  small 
portion  of  the  interaction  region  onto  the  entrance  slit  of  the  vertically 
mounted  SPEX  monochromator,  permitting  a  horizontal  slice  of  the  plasma 
about  5  mm  long  and  50  to  250  pm  wide  to  be  observed.  A  representative 
spectral  scan  of  the  sodium  emission,  extending  from  626  nm  to  452  nm,  is 
presented  as  figure  15.  A  higher  resolution  spectral  scan  of  the  42D-32P 
multiplet  is  displayed  as  figure  16.  In  a  later  experiment  the  focussing 
lens  was  adjusted  so  that  the  region  of  observation  was  scanned  vertically 
allowing  the  lateral  distribution  of  the  emission  to  be  observed. 


Electron  Temperature  Measurements 

The  steady  state  solution  to  the  one  dimensional  radiative  transfer 
equation  enables  the  spectral  radiance  Jnm(v,  *)  arising  from  the  n  to  m 
transition  and  emitted  in  the  x-direction  from  a  slab  of  plasma  to  be 
written  in  the  form 


J  ( V, 
nm' 
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i  x 
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where 
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represents  the  volume  emission  coefficient  with  Nn(x)  the  upper  state 
population  density  at  x,  Anm  the  Einstein  spontaneous  emission  probability 
and  L  (v,x)  is  the  associated  line  profile  function  at  x.  If  the  plasma 
is  assumed  to  be  in  local  thermodynamic  equilibrium  we  may  write 


P(v,  X)  =  [exp{-hv/kTe(x)  -  l}]'1  (33) 


: 


which  is  the  Blackbody  spectral  radiance  (or  Planck  function)  with  hv  as 
the  nm  transition  photon  energy  and  Tg(x)  the  free  electron  temperature  at 
x. 

If  the  plasma  is  uniform  equation  (31)  simplifies  to  the  form 


e  ( v) 
nm'  ' 


Jnm(v’  =  P(v’  Te^1  "  exp^‘  T^.T] 


(34) 


The  optical  thickness  of  the  plasma  at  this  frequency. 
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(35) 
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If  tnm(v»  *)  <<  I*  so  Plasma  is  optically  thin,  then 

J  (v,  Jt)  *  (36) 

nm '  '  nm v  '  '  ' 

or  the  radiance  of  this  line 

•W*)  5  ^  &  NnAn»*  <37> 

In  a  plasma  produced  by  laser  resonance  saturation,  we  can  assume  that 
the  level  populations  are  decribed  by  a  Boltzmann  distribution  down  to  the 
pth_ievel.  Under  these  circumstances  we  can  write 


mV1’ 


v  g  A 

nm3n  nm 


]  -  C 
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e 


(38) 


where  gn  is  the  degeneracy  of  the  nth-level  and  C  is  a  number  that  is  the 
same  for  any  transition  and  is  therefore  not  important.  In  the  optically 
thin  limit  we  can  experimentally  treat  the  multiplet  as  a  single  spectral 
line.  If  ^represents  the  photodetector  signal  arising  from  the  multiplet 
emission,  then  we  can  write 


I  =  R  T  l  J  (1)  (39) 

n  n  L  nm 

where  R  depends  on  the  transition  and  involves  the  spectral  response  of 
n 

the  photodetection  system  vrfiile  T  is  a  constant  of  proportional ity  that  is 
common  to  all  transitions  provided  the  same  photodetector  is  used.  The  sum 
in  (39)  includes  all  the  lines  of  the  multiplet  and  n  is  now  taken  to 
represent  the  multiplet.  Consequently,  if  LTE  is  assumed  to  extend  down  to 
the  resonance  levels  we  can  write 
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Jtn[I  K  3  *  constant  -  E  ,/kT 
L  n  nJ  n2  e 


(40) 


where 

K  =  -(R  v  g  A  )- 1  (41 ) 

n  v  n  nm3n  nm '  '  ' 

Plotting  the  right  hand  side  of  equation  (40)  against  E^  for  the  5s-3p, 
4d-3p,  6s-3p,  5d-3p,  7s-3p  and  6d-3p  multiplets  should  give  a  straight  line 
of  slope  ( - 1  /k T&  ). 

In  figure  17,  the  snail  circles;  (o)  represent  the  results  of  one 
experimental  run  and  the  straight  line  is  derived  on  the  basis  of  a  linear 
regression  analysis.  The  slope  of  this  line  provides  us  with  the  apparent 
temperature ,  Tp  =  4328K.  The  deviation  of  the  experimental  points  from 
this  straight  line  suggests  a  pattern  with  transitions  originating  on  the 
s-levels  lying  above  the  line  and  transitions  originating  on  d-levels  lying 
below  it.  This  scatter  of  the  data  strongly  suggests  that  the  plasma  is 
not  optically  thin  for  spectral  lines  terminating  on  the  32P1/2»  3/2 
resonance  levels. 

If  we  can  no  longer  assume  that  the  plasma  is  optically  thin  then  we 
must  use  equation  (34)  to  define  the  spectral  radiance  at  v  arising  from 
the  nm-transi tion .  In  the  case  of  sodium  the  2S-2P  multiplets  comprise  two 
transitions,  while  the  2D-3P  multiplets  comprise  three.  If  we  assume  that 
the  spectral  lines  are  primarily  electron  impact  Stark  broadened,  then  we 
can  write 
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nm 


(  v-  v  -  S  ) 
v  nm  nm ' 


2 

Y 

'nm 


(42) 


where  y  represents  the  frequency  Stark  IfiHW  and  snm  the  Stark  shift  for 
nm 

the  nm  transition.  Both  are  linearly  proportional  to  the  free  electron 
density. 

It  is  evident  from  equations  (32)  and  (34)  that  evaluation  of  the 
spectral  radiance  requires  knowledge  of  the  upper  state  number  density.  If 
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we  assume  a  Boltzmann  distribution  down  to  the  resonance  level,  laser 
saturation  of  the  resonance  transitions  and  continuity,  we  can  write 


N 


2 


G{N0  -  Ne) 
—1  oG — 


(43) 


where  Nq  is  the  initial  (prior  to  irradiation)  sodium  atom  density,  is 
the  free  electron  density. 


G  =  g2/(92  +  9i)  (44) 

and 

a  =  I  (gi/g2)exP(-Ei2/kTe)  (45) 

represents  the  pseudo-partition  function  —  based  on  the  laser  saturated 
resonance  state  acting  as  a  pseudo-ground  state  divided  by  the  g2. 

The  spectral  radiance  at  v  arising  from  the  transitions  that  comprise 
the  n-mul ti plet , 


hvG(N  -N  ) 

Jn(v.  *)  =  P(v.  Te)[l-exp|-  <#(V>  I»JWv>V«P(-En2^Te)>] 


(46) 


where  the  sum  in  (46)  extends  over  the  transitions  that  constitute  the 

n-multiplet.  It  is  evident  that  Jn(v,  *)  is  3  function  of  three 

independent  variables:  N  ,  T  and  X. 

e  e 

Equation  (46)  was  solved  by  means  of  an  iterative  procedure  that 
converged  on  the  optimum  values  of  (Ng,  Tg  and  1)  for  the  set  of 
experimentally  measured  emission  lines.  The  small  triangles  in  figure  17 
represent  the  computed  values,  derived  in  this  manner  (but  assuming 
aG  «  1),  that  best  fit  the  experimental  data.  The  corresonding  one  slab 
electron  temperature  is  2444K  while  the  degree  of  ionization  is  22.8%. 


This  low  value  for  the  electron  temperature  Is  much  less  than 
predicted  by  our  LIBORS  code.  We  feel  that  part  of  this  difference  may  be 
the  Inadequacy  of  a  one  slab  model  In  allowing  for  optical  depth  effects. 
Indeed,  radiation  trapping  of  -the  resonance  radiation  could  create  a 

|  substantial  halo  of  resonance  state  atoms  that  lie  beyond  the  Ionized 

region  and  absorb  the  observed  spectral  lines. 

In  reality  both  the  electron  density  and  temperature  will  vary 
considerably  across  the  plasma  column  created  by  the  laser  pulse.  The 

Q  radial  distributions  of  both  the  electron  density  and  the  electron 

temperature  can  be  evaluated  by  a  much  more  elaborate  series  of 

measurements.  This  involves  measuring  the  spectral  emission  at  a  number  of 
locations  across  the  radius  of  the  plasma  column. 

^  The  spectral  radiance  at  frequency  v,  arising  from  the  nm-transition 

and  observed  in  the  x-direction  at  some  height  y  above  the  axis  of  a 
cylindrical  plasma  column  of  radius  R  can  be  expressed  in  the  form 


/(R2-y2)  S(R  2-y 2)  e  (v,y,x*)dx* 

J ( V »y )  =  /  „  _.enm(v,y,x)exp{-/  - . VTZ  VT*\ - )dx 
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-✓(R2-y2) 


P(v,y,x*) 


(47) 


See  the  top  left  hand  corner  of  figure  18. 

If  the  plasma  can  be  regarded  as  optically  thin  at  this  frequency  and 
azimuthal ly  symmetric  we  can  spectrally  integrate  over  the  transitions  of 
the  n-multiplet  and  write 

/(R2-y2) 

dn(y)  =  2  /  en(y,  x)dx  (48) 

o 


for  the  radiance  observed  in  the  x-direction  at  height  y  and  arising  from 
the  emission  in  the  n-multiplet,  where 

en(y»x)  =  I  enm(y’x)  (49) 
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The  sum  extending  over  the  transitions  that  comprise  the  n-multiplet. 
Alternatively  we  can  write 
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(y)  =  2  / 
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e  (r)rdr 
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/(r2  -  y2) 


(50) 


by  introducing 
transform  to 
coefficient , 


r2  =  x2  +  y2,  which  can  be  inverted,  using  the  Abel 
provide  the  radial  distribution  of  volume  emission 
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As  indicated  earlier  we  can  relate  Jn(y)  to  the  corresponding  photodetector 
signal  In(y)  through  the  relation 

In(y)  -  RnTJn(y)  (51) 

This  enables  us  to  determine  the  radial  distribution  of  emission  on  the 
n-multiplet  in  terms  of  the  measured  signal  variation  across  the  plasma 
column,  that  is  to  say. 
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If  we  then  compare  the  volume  emission  from  two  multiplets  (n  and  q) 
—  in  practice  the  62D-32P  and  42D-32P  of  sodium,  then  we  can  write 


(53) 
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where  is  a  constant  that  depends  on  the  atomic  parameters  of  the  two 


multiplets.  Consequently,  the  radial  electron  temperature  profile 
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(55) 


can  be  evaluated. 

In  figure  18  the  observed  distribution  of  signals  for  both  the  e^-S2? 
and  42D-32P  multiplets  across  the  plasma  column  are  presented  for  an 
experiment  in  which  the  laser  was  tuned  to  the  32S 1/2-32P 3/2  transition  at 
589.0  nm,  the  incident  energy  was  37  mJ  and  the  peak  sodium  atom  density 
was  estimated  to  be  around  1.5*10 16  cm'3.  The  magnification  of  the  optics 
was  1:1,  with  the  SPEX  monochromator  entrance  slit  set  at  1  cm  *  40  nn. 
The  signal  from  the  RCA  7265  photomultiplier  tube  was  processed  by  the  EGG 
signal  averager  gated  such  that  it  sampled  the  emission  in  a  2  ns  interval, 
approximately  60  ns  after  the  onset  of  the  incident  laser  pulse.  Since  the 
laser  pulse  duration  was  under  40  ns,  the  influence  of  the  laser  field 
(through  the  dynamic  Stark  effect)  should  be  minimal  and  spontaneous  decay 
of  the  resonance  state  population  should  minimize  optical  depth  effects 


(unfortunately  all  of  the  lines  accessible  with  our  present  photomultiplier 
terminate  on  one  of  the  resonance  levels). 

The  curves  drawn  through  the  experimental  points  represent  a  least 
square  fit  using  a  polynomial  of  the  form 

Q(y)  =  {Cx  +  Cjy2  +  C^JexpC-C^y2)  (56) 

where  the  fitting  parameters  Clt  C2»  C4  and  C5  are  chosen  to  minimize  the 
residual  sum  of  the  squares  of  the  deviations  from  the  actual  data  points. 
The  Abel  inversion  of  these  polynomials  is  then  possible  and  the  ratio  of 
the  inverted  curves  is  used  to  estimate  the  radial  electron  temperature 
distribution.  The  experimental  data  presented  in  figure  18  gives  rise  to 
the  radial  electron  temperature  profile  shown  in  figure  19.  Three 
theoretical  temperature  profiles  associated  with  three  very  different  laser 
radial  distributions  are  also  displayed  in  figure  19.  Each  was  synthesised 
from  a  series  of  LIBORS  computer  runs  of  different  incident  laser 
irradiance  as  discussed  earlier.  The  sharp  temperature  peaks  predicted  for 
two  of  the  cases  are  very  unlikely  to  occur  in  reality  due  to  thermal 
conduction.  As  indicated  earlier  we  are  presently  working  on  a  new  code 
that  will  allow  for  thermal  conduction.  Although  this  will  invariably 
smooth  out  the  temperture  profile  -  there  remains  an  important  difference 
in  the  magnitude  of  the  temperature  predicted  by  our  code  and  that 
determined  experimentally.  Currently  we  are  exploring  several  possible 
explanations  for  this  difference. 

It  should  be  noted  that  there  is  reasonable  agreement  between  the 
spatially  averaged  experimentally  measured  Tg  of  figure  19  and  the 
"apparent"  Tg  determined  from  the  Boltzmann  plot  for  the  six  spectral 
lines,  see  figure  17.  As  indicated  earlier  it  is  quite  likely  that  the 
optical  depth  of  the  plasma  could  be  significant  for  transitions 
terminating  on  the  resonance  levels.  This  will  require  correcting  this 
temperature  profile  for  optical  depth  effects. 

To  accomplish  this  we  plan  to  determine  the  radial  electron  density 

profile  by  deconvoluting  the  Stark  broadened  spectra  of  the  42D-32P 

multiplet.  Then  knowing  Nq  and  assuming  LTE  we  can  estimate  the  radial 

resonance  state  distribution  and  correct  both  N  and  T  measurements  for 

e  e 

self-absorption  by  an  iterative  procedure. 


Spectral  scans  of  the  42D-32P  multiplet  were  taken  for  15  positions 
across  the  plasma  column.  The  spectral  resolution  was  estimated  to  be  0.05 
nm  (FWHW)  with  16  laser  shots  being  averaged  to  constitute  one  measurement 
and  128  such  measurements  comprising  one  spectral  scan.  Two  representative 
spectral  scans  of  the  42D-32P  multiplet  are  presented  as  figure  20.  The 
emission  originating  from  the  central  core  of  the  plasma  (y  =  0.25  mm)  is 
seen  to  suffer  a  much  greater  Stark  broadening  (and  red  shift)  than  the 
emission  arising  from  the  more  weakly  ionized  rim  of  the  plasma  column 
(y  =  2.25  mm).  The  laser  in  this  run  was  tuned  to  the  589.0  nm  sodium 
resonance  line  with  an  energy  of  about  37  mJ.  As  before  the  time  of 
observation  was  60  ns  from  the  start  of  the  laser  pulse  and  the  peak  sodium 
atom  density  was  determined  to  be  about  1.5*1016  cm-3. 

At  the  present  time  we  are  developing  a  computer  code  to  perform  the 
deconvolution  of  this  spectral  information.  Generalization  of  equation 
(52)  leads  to 


£„(».  r)  = 


1 

¥7 

n 


R  dl  (v,  y) 


_dy_ 


/(y2-r2) 


(57) 


Since  l_n(v,  r)  is  proportional  to  ep(v»  r)  we  can  use  Stark  broadening  to 
evaluate  Np(r)  given  an  approximate  value  for  Tg(r).  An  optical  depth 
correction  will  be  undertaken  along  the  lines  discussed  earlier  for  the 
case  of  the  electron  temperature  profile. 

Alternatively,  it  is  possible  to  consider  using  emission  on  either  the 
sodium  62P-42S  or  72P-42S  multiplets  as  these  should  not  be  optically 
thick.  To  accomplish  this  we  would  have  to  use  an  RCA  31034 
photomultiplier  which  as  a  GaAs  photocathode  and  a  reasonable  spectral 
response  to  about  850  nm.  Unfortunately,  this  photomultip  lier  has  to  be 
cooled  with  a  carefully  controlled  system  to  maintain  a  constant 
sensitivity.  Currently,  we  are  in  the  process  of  acquiring  such  a  system 
and  plan  to  compare  the  electron  density  and  temperature  profiles  evaluted 
in  this  manner  with  those  determined  using  the  optical  depth  correction 
procedure. 


Laser  Off-Resonance  Maximum  Absorption  (LORMA)-Effect 

During  the  past  year  we  have  undertaken  an  investigation  of  the 

sensitivity  of  laser  ionization*  based  on  resonance  saturation  to  the 

spectral  detuning  of  the  laser  and  have  made  a  surprising  experimental 
observation.  Under  certain  conditions  it  appears  that  maximum  absorption 
of  the  laser  pulse  arises  when  the  laser  is  detuned  by  about  0.4  to  0.5  nm 
to  either  the  blue  of  the  589.0  nm  or  red  of  the  589.6  nm  sodium  resonance 
lines.  This  minimum  in  the  transmitted  energy  has  been  observed  with  both 
the  long  pulse  (1  nsec)  from  our  old  flashlamp  pumped  dye  laser  and  the 

short  pulse  (30  ns)  from  our  new  Nd-YAG  laser  pumped  dye  laser. 

This  laser  off-resonance  maximum  absorption  (LORMA)-effect  is  most 
pronounced  at  higher  densities.  An  example  is  provided  by  figure  21  where 
a  distinct  minimum  in  the  transmitted  laser  energy  is  seen  to  occur  at 
about  0.4  nm  to  the  blue  of  the  589.0  resonance  line  (marked  as  0  on 

wavelength  scale).  The  peak  sodium  atom  density  in  this  run  was  estimated 
to  be  about  4xl016  cm*3.  It  is  also  apparent  that  very  strong  attenuation 
arises  between  the  resonance  lines  and  that  there  might  also  be  a  local 
minimum  at  the  line  centre  wavelength.  However,  more  recent  low  density 
runs  have  not  substantiated  this  latter  observation  and  have  suggested  that 
there  might  even  be  an  asymmetry  in  the  strength  of  the  L0RMA  effect,  with 
more  attenuation  occurring  on  the  red  side  of  the  589.6  nm  resonance  line. 

Simulation  of  this  LORMA-effect  we  felt  would  constitute  a  good  test 
of  our  new  LIB0RS  computer  code.  The  results  of  three  such  detuning 

computer  runs  for  Fp  =  3.75*1025  photons  cm" 2  sec*1  (27  mJ),  2.625*1025  (19 
mJ)  and  2.25*1025  (11  mJ)  are  presented  as  figure  22.  The  peak  sodium  atom 
density  was  taken  to  be  4><1016  cm* 3  to  correspond  to  the  experimental 

conditions  that  applied  to  figure  21.  It  is  immediately  evident  that  a 
well  defined  minimum  in  the  transmitted  laser  energy  is  predicted  to  arise 
for  a  laser  detuning  of  around  0.4  to  0.5  nm.  It  should  be  noted,  however, 
that  no  computer  runs  with  <  0.14  nm  were  undertaken  due  to  the  strong 
laser  coupling  and  the  subsequent  large  increase  in  computational  time 
accompanying  these  small  values  of  laser  detuning.  Nevertheless,  it  is 

obvious  that  our  new  LIB0RS  computer  code  predicts  this  LORMA-effect. 

The  physical  mechanism  responsible  for  this  effect  can  be  understood 
if  we  recall  that  the  dominant  photon  loss  process  is  absorption  by  the 
resonance  transition  and  that  once  the  free  electron  density  exceeds  about 


1015  cm-3  a  dynamic  balance  is  established  between  the  rate  of  laser 
absorption  and  the  rate  of  superelastic  electron  quenching  of  the  laser 
pumped  resonance  state  population,  see  figure  1.  This  suggests  that  the 
rate  of  laser  energy  absorption  will  be  proportional  to  the  product  of  the 
resonance  state  population  and  the  free  electron  density. 

We  propose  that  the  explanation  of  this  LORMA-effect  lies  in  the  fact 
that  a  slightly  detuned  (AX^  *  0.4  nm)  laser  pulse  still  produces  close  to 
full  ionization,  albeit  with  a  much  longer  ionization  time.  As  a 
consequence  the  product  of  the  resonance  state  population  with  the  free 

electron  density  remains  large  for  an  extended  period  of  time  and  this 

maximizes  the  energy  extracted  from  the  laser  pulse  by  electron 
superelastic  collisions.  In  other  words  when  the  laser  is  sharply  tuned  to 
a  resonance  line  the  ionization  proceeds  so  fast  that  the  sodium 

superelastic  electron  laser  energy  converters,  see  figure  1,  cannot  go 
through  many  cycles  before  they  are  destroyed  by  ionization.  On  the  other 

hand  when  the  laser  is  detuned  far  from  the  resonance  line  the  degree  of 

ionization  never  gets  very  large  and  laser  saturation  may  break  down 
reducing,  in  effect,  both  the  free  electron  density  and  the  resonance  state 
population  density. 

It  is  clear  that  there  exists  some  intermediate  detuning  of  the 

laser  for  which  the  number  of  absorption-superelastic  collision  cycles  is  a 
maximum.  This  constitutes  the  detuning  that  extracts  the  maximum  energy 
from  the  laser  pulse.  A  more  detailed  analysis  is  provided  in  Appendix  A. 
We  may  therefore  conclude  that  this  new  observation  represents  further 
evidence  of  the  primary  role  played  by  superelastic  electron  collisions  in 
extracting  energy  from  the  laser  pulse. 

An  interesting  question  is  raised  by  the  LORMA-effect:  if  more  energy 
is  extracted  from  a  laser  pulse  that  is  detuned  by  0.4  nm  compared  to  one 
tuned  close  to  the  resonance  line  -  should  this  additional  energy  not 
appear  in  the  radiation  emitted  by  the  resulting  plasma?  In  figure  23  we 
provide  evidence  of  just  such  an  observation  taken  some  time  ago  with  our 
old  flashlamp  pumped  dye  laser.  The  upper  trace  of  figure  23(a)  displays 

the  incident  laser  pulse  photodiode  signal  while  the  lower  trace  provides 

the  recombination  radiation  in  a  5.0  nm  band  centred  about  405  nm  and 
detected  by  a  photodetection  system  comprising  a  JY-monochromator  and  an 
RCA  1P28  photomultiplier.  In  this  experiment  the  laser  was  tuned  to  the 
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589.6  nm  resonance  line  of  sodium,  while  in  figure  23(b)  the  same 
information  is  displayed  but  for  a  laser  pulse  that  is  detuned  to  590.1 
nm. 

It  is  quite  clear  from  figure  23  that  the  amount  of  radiative 
recombination  emitted  when  the  laser  is  detuned  by  0.5  nm  is  greater  than 
when  it  is  tuned  to  the  589.6  nm  resonance  line.  Close  inspection  of  this 
figure  also  reveals  another  important  difference  that  is  in  keeping  with 
our  theoretical  results  -  there  is  virtually  no  delay  between  the  onset  of 
recombination  radiation  for  the  detuned  laser  while  there  is  a  substantial 
delay  in  the  case  of  the  tuned  laser  pulse. 

Thermal  Conduction  and  LIBORS 

Our  attempts  at  modeling  the  three  dimensional  nature  of  LIBORS  have 
revealed  that  strong  electron  temperature  gradients  would  tend  to  develop, 
see  for  example  figures  9,  10  and  19.  In  reality  the  tendency  to  produce 
strong  electron  temperature  gradients  would  be  countered  by  electron  energy 
transport  through  thermal  conduction.  Our  theoretical  treatment  to  date 
has  omitted  thermal  conduction.  The  justification  for  this  omission  is 
open  to  question  in  light  of  our  recent  findings. 

If  we  take  electron  thermal  conduction  into  account  equation  (14), 
which  describes  conservation  of  electron  energy,  is  modified  and  takes  the 
form 

3u  LIBORS 

Ce  <68> 

where  ue  =  /  7  mgv2f(v)d3v  represents  the  electron  kinetic  energy  density, 

£e  =  /  7  mfiv2vf(v)d3v  represents  the  electron  kinetic  energy  flow  vector 
LIBORS 

and  Cg  stands  for  all  of  the  collisional  energy  sources  and  sinks 

previously  introduced  on  the  right  hand  side  of  equation  (14).  It  should 
be  noted  that  the  form  of  equation  (58)  assumes  that  Joule  heating  is 
negligible.  In  the  absence  of  an  applied  electric  field  this  assumption  is 


reasonable  due  to  ambipolar  effects.  '  The  same  reasoning  permits  us  to 
write, 

9*  =  (59) 

where  X  represents  the  electron  thermal  transport  coefficient. 

If  the  energy  flow  is  not  excessive  the  electron  velocity 
distribution,  f(v),  will  not  deviate  too  greatly  from  the  equilibrium 
electron  velocity  distribution,  f0(v),  and  we  can  express  it  in  the  form 
of  the  truncated  expansion, 

f(v)  =  fQ(v)  +  fj(v)*v/v  (60) 

If  this  is  used  in  the  Boltzmann  equation  for  the  electrons  and  the 
resultant  equation  divided  into  two  differential  equations,  one  of  them 
takes  the  form 
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where  me  is  the  mass  of  the  electron  and  represents  the  electric 
field.  The  right  hand  side  of  equation  (61)  arises  from  a  "Mixture 

Rule", where  represents  an  effective  momentum  exchange  collision 

frequency  that  accounts  for  all  forms  of  elastic  collisions  between  the 
electrons  and  the  heavy  species  such  as  sodium  atoms,  argon  atoms  and 

sodium  ions. 

Now  in  general  the  collision  time  is  considerably  shorter  than  the 

typical  time  for  a  significant  change  in  the  velocity  distribution  and 

consequently  we  may  write 
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If  f o ( v )  is  Maxwellian,  i.e., 
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then  equation  (62)  can  be  expressed  in  the  form 
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If  we  use  equation  (60)  for  f(v)  in  the  definition  of  q  we  arrive  at  the 


expression 
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which  can  be  rewritten  in  the  form 
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where  the  electron  current  density 
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with  AX  =  1.012,  Aa  =  0.476  and  A^1  =  0.  6776.  JtnA  is  the  usual  Coulomb 
logarithm. 

The  first  term  on  the  right  hand  side  of  equation  (76)  represents  the 
sum  of  the  electron  momentum  transfer  collision  frequencies  for  each  of  the 
sodium  atomic  states  (primarily  ground  and  resonance  state  atoms)  and  the 
argon  atoms,  while  the  second  term  represents  the  electron-ion  momentum 
transfer  collision  frequency.  The  three  different  coefficients  AX,  Act  and 
A1*1  are  required  in  order  that  the  expressions  used  agree  in  both  the  low 
and  high  ionization  limit  with  those  obtained  when  electron-electron 
collisions  are  rigourously  included. 

If  the  plasma  is  fully  ionized  then  the  first  term  on  the  right  hand 
side  of  equation  (76),  representing  the  sum  of  collision  frequencies 
involving  neutral  atoms,  is  zero  and  X  becomes  the  Spitzer-Harm  value,  X^. 
Ambipolar  effects  are  found  to  make  the  first  term  on  the  right  hand  side 
of  (70)  negligible  compared  to  the  second,  so  that  we  usually  write 

S*  -  -*VTe  (77) 

It  should  be  noted  that  X  is  in  general  smaller  than  X^H  and  it  therefore 
follows  that  electron  thermal  conduction  for  a  partially  ionized  vapor  is 
less  than  for  a  fully  ionized  vapor. 

For  the  case  of  a  sodium  plasma  formed  by  laser  resonance  saturation 
the  ratio  of  X/X<-H  has  been  computed  for  an  original  sodium  density  (prior 
to  laser  excitation)  of  1016  cm" 3  as  a  function  of  the  free  electron 
temperature.  The  results  are  presented  in  figure  24. 

The  full  curve  ( - )  corresponds  to  the  case  where  LTE  is  assumed  down 

to  the  resonance  state  as  the  resonance  transition  is  laser  saturated, 

while  the  broken  curve  ( - )  corresponds  to  LTE  down  to  the  ground 

state.  Also  displayed  by  the  dotted  curve  ( . )  in  figure  24  is  the 


variation  of  with  Tg  for  the  situation  where  argon  at  a  density  of 
5*1016  cm  ‘3  is  added  to  the  first  case.  The  variation  of  with 
electron  temperature  is  also  provided  by  the  ( — A — )  curve  in  figure  24. 

These  results  clearly  indicate  that  the  Spitzer-Harm  electron  thermal 
conductivity  is  reliable  to  within  15%  of  the  exact  value  for  electron 
temperatures  greater  than  4000  K,  irrespective  of  the  model  used.  For 
electron  temperatures  less  than  this  value  the  thermal  conductivity  rapidly 
drops  and  we  could  probably  ignore  electron  thermal  conduction. 


Self-Pumped  Resonance  Inversion  Scheme 

An  intense  flash  of  light  was  used  to  excite  the  first  laser  and  many 

of  today's  most  popular  lasers  involve  optical  pumping.  In  the  race  to 

f  29-34 1 

develop  an  X-ray  laser,  radiative  pumping  schemes  figure  prominently:  ‘ 

In  almost  all  of  these  a  coincidence  is  sought  between  the  wavelength 

emitted  by  the  radiating  ion  species  and  that  absorbed  by  the  recipient  ion 
(34) 

species.  ' 

In  this  section  we  consider  a  new  approach  towards  the  development  of 
short  wavelength  (EUV  to  X-ray)  lasers  based  on  sudden  radiative  pumping. 
The  central  feature  of  this  "self-pumped  resonance  inversion"  (SPR I  )-scheme 
is  the  large  fine  structure  splitting  of  the  ground  state  found  in  members 
of  the  thallium  (and  to  a  lesser  extent  in  the  indium-  )-i  soelectronic 
series. 

In  order  to  explain  the  SPRI-scheme  we  shall  consider  a  representative 
example  of  the  TX-i soel ectronic  series. . .doubly  ionized  bismuth.  A  partial 
Grotrian  energy  level  diagram  of  the  lowest  BiHI  levels  is  presented  as 
figure  25.  If  a  plasma  containing  Billl  as  a  major  constituent  is  suddenly 
heated  a  burst  of  emission  on  all  the  strong  lines  will  ensue.  If  this 
radiation  is  directed  at  a  cold  Billl  plasma  photoexcitation  from  the 
6 2 P 1/2  ground  level  will  rapidly  populate  many  of  the  higher  levels, 
including  the  72Sj/2  and  62D3/2  levels.  The  upper  level  of  the  split 
ground  state,  62P3/2,  will  not,  however,  be  excited  in  this  manner  and  this 
suggests  that  it  may  be  possible  to  create  a  transient  population  inversion 
between  the  62D3/2  (or  to  a  lesser  extent  the  72S3/2)  and  62P3/2  levels. 
This  concept  is  schematically  illustrated  in  figure  26. 

Clearly,  the  larger  the  population  initially  residing  in  the  62P3/2 
level  the  more  difficult  it  will  be  to  create  the  population  inversion. 
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This  tends  to  impose  a  restriction  on  the  cold  plasma  temperature  for  a 
given  source  radiance  or,  equivalently,  the  minimum  jump  in  the  temperature 
of  the  radiating  plasma. 

Although  there  are  many  ways  to  rapidly  heat  a  plasma,  the  temperature 

achieved  is  usually  not  controlled  and  uniformity  rarely  attained.  Laser 

saturation  of  a  low  lying  resonance  transition  of  an  ion  that  forms  a  major 

(351 

constituent  of  a  plasma  does,  however,  offer  just  these  features'  '  by  way 
of  superelastic  electron  heating. 

The  controlled  jump  in  the  electron  temperature  predicted  for  this 
superelastic  heating  technique  is  a  direct  result  of  the  self-limiting 
nature  of  the  interaction.  In  essence,  once  laser  saturation  of  the 
resonance  transition  is  achieved,  the  electron  temperature  is  determined  by 
a  dynamic  balance  of  superelastic  and  inelastic  collisional  processes  and 
is  therefore  independent  of  the  laser  irradiance  or  even  the  density  of  the 
relevant  ion.  Such  a  controlled  jump  in  the  free  electron  temperature  is 
important  from  both  the  standpoints  of  efficiency  and  optimization  of  the 
spectral  composition  of  the  emitted  radiation. 

In  general,  the  radiating  ion  used  in  the  SPRI-laser  scheme  will  not 
be  suitable  for  superelastic  electron  heating  and  as  a  consequence  the  hot 

plasma  emission  source  will  probably  have  to  comprise  both  kinds  of  ions: 

those  involved  in  heating  the  free  electrons  and  those  that  emit  the 

appropriate  line  radiation  for  pumping  the  cold  plasma. 

In  our  preliminary  discussion  of  the  self-pumping  resonance  inversion 
scheme  we  used  the  model  indicated  in  figure  26.  In  essence  a  plasma 

comprising  doubly  ionized  bismuth  (Billl)  was  assumed  to  be  suddenly  heated 
such  that  the  free  electron  temperature  jumped  to  a  value  Tg  in  a  time  that 
is  short  compared  to  times  of  interest.  A  subsequent  step-like  pulse  of 
104  nm  line  radiation,  corresponding  to  the  (62D3/2  -  6  2P  i  / 2 )  transition 
(see  figure  25)  would  be  emitted.  This  flash  of  radiation  would  propagate 
into  the  cold  low  density  plasma  of  Billl  and  serves  to  excite  the  62D3/2 
level  of  the  recipient  ion.  Since  there  would  be  no  corresponding 
excitation  of  the  upper  level  of  the  split  ground  state,  62P3/2,  we  showed 
that  a  population  inversion  was  created  between  levels  62D3/2  and  62P3/2, 
given  a  small  enough  initial  population  on  the  62P3/2  level. 

In  this  analysis  we  neglected  the  (72S  -  62P)  multiplet  (see  figure 


25)  as  these  transitions  are  somewhat  weaker  and  we  also  neglected  the 
effect  of  radiation  at  122.5  nm  correspondi ng  to  the  (62P3/2  -  6%5/2) 
transition.  We  further  assumed  that:  the  density  in  the  cold  plasma  was 
low  enough  that  collisional  mixing  and  broadening  could  be  neglected,  while 
that  of  the  hot  plasma  was  adequate  to  ensure  near  blackbody  emission  on 
the  lines  of  interest. 

Under  these  circumstances  we  were  able  to  show  that  a  population 
inversion  could  be  created  for  a  hot  plasma  temperature  of  9  eV  and  a  cold 
plasma  temperature  of  1  eV.  In  the  case  of  a  cold  plasma  density  of  1014 
cm"3  a  peak  inversion  density  of  several  percent  was  predicted  to  arise  on 
the  132.7  nm  transition  within  a  distance  of  2*10"2  cm  of  the  sharp 
boundary,  see  figure  27.  This  small  gain  region  is  a  direct  result  of  the 
short  optical  path  of  the  104  nm  pump  radiation  and  this  has  prompted  us  to 
consider  if  it  would  not  be  more  fruitful  to  consider  the  population 
inversion  that  would  be  created  on  the  134.6  nm  (72S1/2  -  62P3/2)  line 
through  pumping  on  the  weaker  105.2  nm  line. 

Within  the  past  few  months  we  have  undertaken  a  more  thorough  analysis 
of  this  concept.  This  treatment  takes  account  of  all  of  the  radiative 
processes  coupling  the  72S,  62D  and  62P  terms  of  BiHI  and  involves  solving 
the  following  set  of  population  rate  equations: 
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for  the  five  levels  indicated  in  figure  25. 

If  we  neglect  higher  states  of  excitation  then  continuity  can  be 
expressed  by  the  relation 


where  N0  represents  the  initial  Billl  ion  ground  state  density,  prior  to 
irradiation.  A  further  simplification  is  afforded  by  assuming  that  the 
population  of  levels  4  and  5  are  always  in  the  ratio  of  their  degeneracies 
—  this  should  be  a  conservative  assumption.  These  two  assumptions  enable 
the  set  of  five  differential  equations  indicated  by  equation  (78)  to  be 
reduced  to  three  population  equations. 

If  the  cold  Billl  plasma  is  assumed  to  be  suddenly  irradiated  by 

U 

blackbody  radiation  at  temperature  T  ,  then  the  coefficients  in  the 
population  rate  equations  are  constants  and  solutions  can  be  obtained 
directly.  Since  we  are  assuming  that  the  density  of  the  plasma  is 
sufficiently  low  (1014  cm-3)  that  collisional  effects  are  negligible  on  the 
time  scales  of  interest  —  we  shall  consider  the  lines  to  be  Doppler 
broadened. 

This  more  complete  analysis  reveals  that  no  population  inversion  is 
created  on  either  the  72S i/2  -  62P3/2  or  672D3/2  -  62P3/2  transitions  for 
the  hot  and  cold  plasma  temperatures  of  9  eV  and  1  eV,  respectively.  It 
should  be  noted  that  the  temperature  of  the  cold  plasma  only  enters  in 
determining  the  initial  population  of  the  62P3/2  terminal  level  and  the 
1  eV  was  felt  to  be  the  minimum  that  would  make  sense  for  a  Billl  plasma 
that  was  created  sometime  before  irradiation.  If,  however,  the  Billl 
plasma  was  generated  very  rapidly  just  prior  to  irradiation  (maybe  through 
multiphoton  ionization  )  it  may  be  possible  to  have  a  low  62P3/2  level 
population  —  equivalent  to  a  0.25  eV  Boltzmann  distribution.  We  have 
found  that  under  these  circumstances  a  population  inversion  of  0.5%  and 
1.2%  can  be  created  on  the  62D3/2  "  62f>3/2  and  ^2Si/2  “  ^ 2p  3  /  2 
transitions,  respectively.  It  should  also  be  noted  that  our  assumption 
concerning  the  coupling  of  the  populations  on  the  62D5/2  and  62D3/2  levels 
may  be  too  conservative  since  the  collisional  relaxation  time  is  likely  to 
be  much  greater  than  the  times  of  interest.  If  we  decoupled  these  two 
levels  we  would  probably  obtain  inversions  for  somewhat  higher  effective 
temperatures  of  the  cold  plasma. 

New  Three  Photon  Saturation  Laser  Scheme 

During  the  past  six  months  we  have  been  studying  the  feasibility  of  a 
new  approach  towards  the  development  of  a  short  wavelength  laser  that  is 
based  on  three  photon  saturation  of  an  ion  possessing  a  low  lying 
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metastable  state.  The  basic  inversion  scheme  for  the  case  of  singly 
ionized  strontium  Sr  1 1  is  indicated  in  figure  28.  In  essence  we  propose 
that  a  momentary  population  inversion  should  be  created  between  the  52F  and 
42D  terms  if  the  former  is  directly  populated  from  the  ground  state  by  a 
three  photon  process  while  the  latter  is  minimally  excited. 

The  key  to  this  concept  is  the  rapid  formation  of  the  Sr  1 1  ground 
state  ions  through  two  photon  ionization  of  strontium  vapor,  and  the 
simultaneous  three  photon  excitation  of  the  52F  term  by  the  same  laser 
field.  This  should  provide  a  means  of  populating  the  52F  term  before  any 
appreciable  population  is  produced  on  the  42D  metastable  term.  The 
wavelength  of  the  radiation  (422  nm)  required  for  this  scheme  is  close  to 
the  resonance  radiation  of  the  Sr  1 1  ion  and  as  a  consequence  the  three 
photon  excitation  cross-section  should  be  considerably  enhanced  by  this 
near  resonance  condition.  Furthermore,  the  presence  of  the  52D  term  should 
also  improve  the  strength  of  this  interaction.  If  this  concept  is 
successful  it  could  lead  to  an  effective  method  of  obtaining  laser 
radiation  at  close  to  the  third  harmonic  of  the  pump  laser  wavelength  — 
but  with  much  higher  efficiency  than  can  be  achieved  through  present 
methods  of  third  harmonic  generation.  It  could  also  be  extended  to  other 
wavelengths  by  using  different  transitions  and  other  ions  possessing  low 
lying  metastable  states. 

In  order  to  test  the  feasibility  of  this  method  of  creating  a  short 
wavelength  laser  we  have  initiated  a  combined  theoretical  and  experimental 
program.  These  programs  are  designed  to  answer  three  important  questions. 
First  can  we  two  photon  ionize  a  reasonable  length  of  strontium  vapor 
with  our  Nd-YAG  laser  pumped  dye  laser,  secondly  can  we  excite  a  sufficient 
number  of  the  Sr  1 1  ground  state  ions  into  the  52F  term  to  create  a 
population  inversion  between  52F  and  42D  terms  and  thirdly  is  competition 
from  either  of  the  other  52F-n2D  (n  =  6  or  5)  multiplets  likely  to  prevent 
laser  action  on  the  52F-42D  multiplet  at  177  nm. 

Theoretical  Program 

We  have  modeled  the  Sr  1 1  ion  by  the  5-level  system  shown  as  figure  29. 
It  is  evident  that  laser  radiation  at  about  422  nm  can  populate  level  5 
(corresponding  to  the  52F  term  of  Sr II)  through  three  photon  near  resonant 
absorption.  If  the  ions  are  created  cold  in  the  ground  state  at  virtually 
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the  same  moment  via  two  photon  ionization  of  strontium  vapor  by  the  same 
laser  pulse  then  clearly  a  population  inversion  would  be  created  on  the  5-4 
and  5-3  transitions  in  addition  to  the  desired  inversion  on  the  5-2 
transition. 

The  rate  of  three  photon  excitation  can  according  to  Yariv^6^  be 
written  in  the  form 


WP 


27ne6IL3 


g15flrfr3f35 


32(meeocfi)3  <^3^3 >  2^31~Zui )  2au,51 


(80) 


where  f  represents  the  absorption  oscillator  strength  of  the  n-m 
mn  r  r  3 

transition  and  u>nm  the  corresponding  angular  frequency.  The  resonance  and 
intermediate  levels  involved  in  the  three  photon  process  are  designated  "r" 
and  "3"  respectively.  Au>51  represents  the  bandwidth  and  g15  the  degeneracy 
ratio  (Qi/gs)  of  the  5-1  transition.  I[_  is  the  laser  irradiance  and  the 
other  constants  have  their  usual  values. 

In  reality  there  are  several  three  photon  paths  that  would  lead  to 
excitation  of  the  52F  term: 

(i)  52S1/2  52P1/2  -*•  52D3/2  82f:5/2 

•  (31)  82S1/2  *  ^  2^  3  /  2  +  ^ 3  /  2  *  ^2F5/2 

(iii)  52S1/2  52P  3/2  52°5/2  *  5*F5/2 

(iv)  52S1/2  *  8  2P  3  /  2  +  5  2°  5/2  52F7/2 

If  the  laser  wavelength  is  tuned  such  that  three  photons  just  excite  the 
52F  term  (i.e.,  X^  =  422.14  nm)then  the  near  resonance  with  the 
52S i / 2~5 2P 1/2  transition  (at  421.67  nm)  ensures  that  we  need  only  consider 
the  first  path  as  its  rate  will  exceed  the  rates  associated  with  the  other 
paths  by  nearly  two  orders  of  magnitude  due  to  the  small  value  of  the 
(w  .-ai)  factor  in  equation  (80). 


Under  these  circumstances  the  three  photon  excitation  rate  per  ion 


I.  3(t) 

Wn(t)  *  1. 96x10“ 18  -i - 

P'  '  Au51 


where  I L  ( t )  is  the  laser  irradiance  at  time  t  and  Aw51  is  the  bandwidth  of 
the  three  photon  transition.  The  appropriate  population  rate  equations  for 
our  5-level  model  are: 


-rr-  =  N.W  -  Nr.Rr,  -  N„R„  -  N„Rr,  -  N(.Ar 
dt  1  p  54  54  53  53  52  52  5  5 


-rr—  =  Nc  Rc  +  N,AC  -  N  A 
dt  5n  5n  5  5n  n  n 


where  n  =  4,  3  and  2.  In  general,  is  the  population  density  of  the  m 

level,  and  N  =  N  -  g  N  represents  the  population  inversion  density 
nm  n  Jnm  m 

between  levels  n  and  m.  The  stimulated  emission  rate 


R5n  5  7T1  J  I5n(v)L5n(v)dv 


where  B^n  represents  the  Milne  stimulated  emission  coefficient,  ( v)  the 
line  profile  function  and  I5n(v)  the  spectral  irradiance  that  arises  from 
the  5-n  transition.  The  total  rate  of  spontaneous  decay  of  level  n. 


where  the  sum  extends  over  all  levels  m  lower  than  n. 

If  we  assume  that  the  radiation  field  originating  from  stimulated 
emission  on  the  5-n  transition  (where  n  =  4,  3  and  2)  is  much  narrower  than 
the  appropriate  atomic  linewidth  and  centred  at  the  appropriate  line  centre 
frequency,  v^,  then  we  can  write 

15„<V)  -  I5n6(v  -  v5n)  (86) 

where  I5n  is  the  irradiance'  of  the  radiation  field  and  6(  v  -  v^n)  the 
appropriate  delta  function.  Under  these  conditions 


r  =  i  I 

K5n  47:  A5nL5n 


where  Lc  represents  the  value  of  the  line  profile  function  at  the  line 
bn 

centre  frequency  v^. 

In  order  to  calculate  the  irradiance  produced  on  each  of  the  inverted 
transitions  along  the  path  of  the  pump  beam  we  must  obviously  solve  the 
above  set  of  rate  equations  combined  with  the  appropriate  radiative 
transfer  equations.  In  the  first  place  we  must  take  account  of  the 
attenuation  of  the  pump  laser  pulse.  Since  the  speed  of  light  can  be 
regarded  as  effectively  infinite  for  our  conditions  and  the  laser  beam  is 
assumed  to  propagate  primarily  in  the  z-direction  we  can  write 


ai.  (z,t) 

- ^ -  =  -N1(z,t)piL3(z,t)  -  |Nq  -  N1(z,t)}aIL2(z,t) 


where  P  represents  the  three  photon  excitation  rate  l**^*1^ - }  *  3hvL,  a 

the  two  photon  ionization  rate  *  2hv[_,  N0  the  strontium  atom  density 
prior  to  laser  irradiation,  Ni(z,t)  the  strontium  ion  ground  state  density, 
and  Ii(z,t)  the  pump  laser  irradiance.  hvL  represents  the  laser  photon 


energy.  Equation  (88)  has  to  be  solved  in  conjunction  with  the  ion 
creation  equation: 


dN^z.t) 
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(89) 


We  shall  now  consider  the  growth  of  the  radiation  field  on  each  of  the 
5-4,  5-3  and  5-2  transitions.  To  describe  the  amplification  of  the 

radiation  field  on  the  5-n  transition  propagating  in  the  z-direction  we 
again  use  the  steady  state  one  dimensional  radiative  transfer  equation  — 
which  in  this  instance  takes  the  form 


9_ 
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(v.z)} 


hv5nA5nN5L5n(v)  Tn  +  °5n ( v^N5nI5n ( v,z ) 


(90) 


where  I5n(v,z)  represents  the  spectral  irradiance  at  z  arising  from  the  5-n 

h  vs 

transition  and  o5n(v)  I2  B5nL5n(v^  represents  the  corresponding 

stimulated  emission  cross  section  with  n  =  4,  3  and  2. 

If  we  integrate  over  a  range  of  frequencies  that  sufficiently 
encompass  the  5-n  line  and  takes  account  of  our  assumption  regarding  the 
narrow  bandwidth  of  the  emission  we  can  write 


dI5nU) 


Hz-  =  E5nAQ  +  “Sn^n^ 


(91) 


n5a 

emission  coefficient,  lgn(z)  is  the  irradiance  of  the  radiation  at  position 
z  and 


}  is  the  spectrally  integrated  volume  spontaneous 


where  e 
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(92) 


is  the  gain  coefficient  at  the  5-n  line  centre  frequency,  v^n.  As  before 
Ngn  represents  the  population  inversion  density  on  the  5-n  transition. 

If  we  assume  that  the  gain  medium  (the  strontium  plasma)  is  reasonably 
uniform  over  the  spatial  range  of  interest,  then  the  solution  of  equation 
(91)  is  straightforward  and  we  can  write 

fltr  Z  £r  AQ  Ctr  Z 

i6n(z>  ■  +  -42—  <M> 


where  1^(0)  represents  the  irradiance  (arising  from  emission  on  the  5-n 
transition)  incident  on  the  segment  of  plasma  under  consideration.  Clearly 
enumeration  of  equation  (93)  is  required  in  order  to  provide  the 
appropriate  stimulated  rates,  equation  (87),  that  arise  in  the  population 
rate  equations.  However,  the  equations  are  coupled  as  e^n  and  depend 
upon  the  appropriate  population  densities. 

It  is  also  apparent  from  equations  (87)  and  (92)  that  we  need  to 


evaluate  the  effective  bandwidth  of  the  5-n  transition,  Ac  (=  1/L*- _ }  * 

f371  3n  3n 

which  in  turn  depends  upon  the  gain  of  medium.  Yariv'  '  indicates  that  in 


the  case  of  an  inhomogeneously  broadened  line  the  effective  bandwidth  for 
the  unsaturated  regime  (i.e.,  I^n  <  I^n}  is  given  by 


where  Avl?  is  the  FWHM  Doppler  width  for  the  5-n  transition  and  icc  is  the 
5n  bn 

incremental  line  centre  gain  constant  in  the  limit  of  extreme  inhomogeneous 
broadening  and  is  defined  by  the  relation 
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(95) 


N5°n 

4n3/2v§n  Av5n 


Here  N^n  {=  W p Nq / A&  }  represents  the  unsaturated  population  inversion 
density  for  the  5-n  transition.  The  saturated  irradiance 


°5n  T5n 


where  a,  is  the  5-n  line  centre  stimulated  emission  cross-section  and 

n 

T5n  5  ly/A5n"  In  the  saturated  regime  { i .  e . ,  I^n  >  I^n}  Yariv^37^  shows 

that  A,  =  Avc  . 

5n  5n 

In  a  first  attempt  to  model  the  growth  of  the  radiation  field  on  each 
of  the  5-n  transitions  (n  =  4,  3  and  2)  resulting  from  three  photon 

excitation  of  level  5,  we  have  assumed  that  two  photon  ionization  of 
strontium  vapor  produces  a  uniform  column  of  ground  state  SrII  ions.  We 
further  assumed  that  negligible  laser  power  was  absorbed  in  this  process. 
The  justification  for  this  lies  in  the  small  amount  of  laser  energy 
required  for  this  purpose  compared  to  that  available.  Although  this 
calculation  was  somewhat  simplified  it  shows  that  this  assumption  is 

reasonable  for  our  preliminary  model  of  the  interaction. 

We  then  divide  this  plasma  column  into  a  number  of  thin  slabs  and 

assume  that  within  each  slab  conditions  are  uniform  and  equations  (82)  and 
(83)  can  be  solved  to  describe  the  temporal  variation  of  the  level 

population  densities  subject  to  the  radiation  field  present.  It  should  be 
noted  that  we  have  implicitly  assumed  that  the  plasma  density  is 

sufficiently  low  that  we  can  neglect  collisional  processes  for  the  times  of 
interest.  Of  particular  concern  is  electron  excitation  of  the  4^ 
metastable  term  of  SrII.  In  the  case  of  a  30  ns  laser  pulse  this  requires 

the  electron  density  to  be  no  greater  than  1015  cm-3.  If  full  ionization 

were  achieved  this  implies  that  the  strontium  atom  density  prior  to 

Irradiation  should  be  limited  to  about  1015  cm-3. 
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The  irradiance  arising  from  the  5-n  transition  and  emitted  in  the 
z-direction  from  the  ith  slab  of  thickness  A 2  is  given  by 


lln  e 


“Sn*2 


+  [e«4n* 


where  I5”*  represents  the  incident  irradiance,  e^n  the  volume  emission 
coefficient  and  al  the  gain  coefficient  (corresponding  to  the  5-n 
transition)  for  the  i*"  slab.  In  the  case  of  the  first  slab  (i=l)  we  have 


1  1  A 
E5nAQ  r  “bn*2 

— 1 —  le 


where  AQ  was  set  equal  to  the  area  of  the  ioni2ed  column  divided  by  the 
square  of  the  column  length. 

The  laser  irradiance  used  for  our  preliminary  calculation  approximated 
the  experimental  laser  pulse  and  was  given  by  the  following  expression: 


t-Xj  2 


T  f  f  )  -  T  ^  ^  . 

(t)  *1  6 


where  I  L  represents  the  peak  irradiance,  xx  =  10  ns  and 
_  —  f  6  n  S  f 0  r  t  ^  10  ns  A, fAnn.fiAn  Af  f  Ua  I  . aa a  paHi  af^AA  a  p i r i n  n  f r Am 


42  ns  for  t  >  10  ns 


.  Attenuation  of  the  laser  radiation  arising  from 


the  near  resonant  three  photon  pumping  of  the  strontium  ion  was  also  taken 
into  account  through  the  relation 


If  we  assume  N1  and  g  are  uniform  over  the  i 
solution  of  this  equation  yields 


slab  thickness  Az,  the 


th 
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(101) 


The  above  set  of  equations  has  been  solved  by  a  4th-order  Runge-Kutta 
routine  to  yield  the  emission  on  each  of  the  52F-62D,  5^-5 20  and  52F-42D 
multiplets.  The  results  of  our  preliminary  calculations  indicate  that  the 
equations  are  stiff  and  that  a  better  computational  routine  DEGEER  should 
be  employed.  Nevertheless,  our  initial  calculations,  performed  for  a 
column  length  of  0.2  rmi,  indicated  appreciable  emission  on  the  2.82  pm  and 
177  nm  transitions  as  seen  from  figure  30.  These  results  look  encouraging 
for  they  suggest  that  laser  action  on  the  VUV  line  should  be  attainable. 
Currently  we  are  improving  our  computer  code  in  order  to  extend  these 
conuptations  to  more  practical  lengths. 

Experimental  Program 

In  parallel  with  this  theoretical  work  we  have  designed  and  assembled 
an  experimental  facility  to  both  test  this  concept  and  study  laser 
superelastic  heating  of  an  ionic  resonance  transition  within  a  plasma.  The 
former  will  be  achieved  by  tuning  the  laser  to  422.14  nm,  while  the  latter 
will  be  accomplished  by  tuning  the  laser  close  to  the  421.67  nm  resonance 
line  of  Sr  1 1 .  In  both  cases  the  side  light  emission  will  be  studied,  but 
for  the  three  photon  experiment  we  also  plan  to  monitor  the  emission 
coaxial  with  the  pump  laser  beam. 

The  basic  components  of  this  facility  are  shown  in  figure  31  and 
include  a  crossed  heat  pipe,  with  four  optical  ports,  that  is  designed  to 
operate  with  strontium  vapor.  A  vacuum  blind  phototube  with  a  visible 
wavelength  blocking  filter  will  be  used  to  monitor  the  vacuum  ultraviolet 
emission  at  177  nm  while  a  photomultiplier  coupled  to  a  SPEX  17001 1 


monochromator  will  be  used  to  monitor  the  yellow  562.4  nm  emission.  Our 
new  Nd-YAG  laser  pumped  dye  laser,  will  be  used  to  pump  the  strontium  vapor 
for  both  experiments.  The  signals  from  our  photodetectors  will  be 
processed  by  our  gated  EGG  signal  averager. 

At  the  time  of  writing  this  facility  has  just  become  operational  and 
we  have  undertaken  a  preliminary  experiment  that  involved  scanning  the  dye 
laser  from  0.3  nm  to  the  red  of  422.14  nm  through  to  0.3  nm  to  the  blue  of 
421.67  nm.  The  axial  emission  at  562.4  nm  was  monitored  during  this 
experiment.  Although  a  reasonable  signal  was  observed  it  appeared  that  the 
vapor  density  was  too  high  and  the  three  photon  generated  emission  at  562.4 
nm  was  mixed  into  the  tail  of  what  appeared  to  be  superel asti cal ly 
generated  plasma  emission  associated  with  laser  saturation  of  the  421.67  nm 
resonance  line.  These  results  are  nevertheless  encouraging  for  they 
suggest  that  both  phenomena  could  indeed  be  studied.  We  are  next  planning 
a  lower  density  run  to  optimize  the  conditions  for  the  3  photon  scheme. 
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NEW  DISCOVERIES  STEMMING  FROM  RESEARCH 


Saturation  of  an  atomic  transition  by  the  intense  radiation  field  of 

a  suitably  tuned  laser  represents  an  important  kind  of  interaction  with  a 

wide  range  of  potential  applications.  The  consequences  of  laser  resonance 

saturation  and  the  applications  stemming  from  this  interaction  depend  to  a 

very  large  extent  upon  the  period  of  saturation.  If  the  resonance  to 

ground  level  populations  are  only  momentarily  locked  in  the  ratio  of 

respective  degeneracies  (saturation)  the  principal  effect  is  a  burst  of 

intensified  spontaneous  emission  that  can  be  used  to  diagnose  the  excited 
(1  2) 

medium:  ’  On  the  other  hand  an  extended  period  of  saturation  (lasting 

for  much  longer  than  the  resonance  state  lifetime)  can  lead  to  extensive 

(3  4) 

perturbation  of  the  medium:  '  Indeed,  if  the  free  electron  superelastic 

collision  time  is  short  compared  to  the  duration  of  the  laser  pulse  near 

(4-191 

total  ionization  of  the  laser  pumped  species  is  rapidly  achieved:  1  In 

the  case  of  an  ionic  species  extremely  rapid  changes  of  electron  temperature 
can  be  produced  by  this  means  f2^ 

The  author  was  the  first  to  recognize  the  importance  of  laser  satura- 

(1-41 

tion  and  many  of  its  possible  areas  of  application.  '  Subsequent  work 
in  the  author's  laboratory  and  elsewhere  have  proven  that  laser  resonance 
saturation  does  represent  a  significant  form  of  interaction  between  laser 
radiation  and  atomic  vapours  or  plasmas. 

Momentary  laser  saturation  represents  a  powerful  diagnostic  technique 

that  is  finding  application  in  many  areas  ranging  from  fusion  reactor 

(211  (22  231 

studies  J  to  combustion  measurements:  ’  *  The  combination  of  laser 

ablation  and  laser  saturation  spectroscopy  represents  a  new  approach  at 

evaluating  fundamental  atomic  quantities  such  as:  radiative  lifetimes, 

branching  ratios,  transition  probabilities  and  selected  collision  cross- 

sections.  A  preliminary  paper  on  this  subject  was  published  by  us  in 
(24) 

Physical  Review:  This  technique,  as  well  as  being  convenient  and 

accurate,  is  particularly  well  suited  for  measurements  on  short  lived, 
highly  ionized  species  created  by  laser  ablation.  Furthermore,  it  is 
versatile  and  can  use  multiphoton  or  stepwise  excitation  as  the  means  of 
generating  the  bursts  of  intensified  emission. 

As  a  spin-off  of  this  work  we  have  also  shown  that  this  concept  can 
also  form  the  basis  of  a  new  form  of  trace  element  laser  microprobe  called 
a  TABLASER?255 


More  recently,  we  have  demonstrated  in  a  proof  of  principle  experiment 

that  laser  saturation  spectroscopy  can  be  used  to  directly  measure,  with 

both  spatial  and  temporal  resolution,  the  ion  to  neutral  atom  density  ratio 

in  a  rapidly  expanding  ablation  plasma .  ^ 

The  ionization  capabilities  of  extended  laser  resonance  saturation 

have  now  been  indisputably  demonstrated  with  a  variety  of  experiments^ 

We  have  developed  a  model  of  this  laser  ionization  based  on  resonance 

saturation  (LIBORS)  and  have  shown  that  the  interaction  can  be  thought 

to  proceed  in  four  stages.  ’  *  Confirmation  of  certain  aspects  of  our 

theory  has  been  provided  by  several  research  groups^  Recently,  there 

has  been  some  success  in  using  laser  resonance  saturation  for  the  purpose 

f  271 

of  creating  infrared  laser  action.  '  Currently,  our  efforts  are  concen¬ 
trating  upon  obtaining  a  better  understanding  of  this  interaction  and  in 
studying  how  best  to  employ  this  fast  and  very  efficient  method  of  coupling 
laser  energy  into  a  plasma  for  the  purpose  of  developing  an  efficient  short 
wavelength  laser. 
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FIG.  1  SUPERELASTIC  ELECTRON  HEATING  THROUGH  LASER  RESONANCE 
SATURATION  OF  SODIUM  ATOMS. 


FIG.  2  BASIC  COLLISIONAL-RADIATIVE  MOOEL  USED  IN  5  LEVEL-LIBORS 
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FIG.  3  COMPARISON  OF  20-  AND  5-lEVEt  LIBORS  COMPUTER  CODES 
FOR  SODIUM  ATOM  DENSITY  OF  10u  CM'*. 

(•)  Temporal  Variation  of  N2/N.  for  the  Truncated  Laser 
Pulse  Indicated  by  F  (Photons  cm'2  sec'1)* 

(b)  Temporal  Variation  of  Ne,  Te  and  Tj  for  the  same  F. 

(c)  Temporal  Variation  of  Nj^  for  the  same  F. 
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FIG-  4  COMPARISON  OF  20  AND  5-LEVEL  LIBORS  COMPUTER  CODE 
PREDICTIONS  OF  TEMPORAL  VARIATIONS  IN  Ne,  Te,  T1  and 
N^/N^  FOR  LASER  PHOTON  FLUX  DENSITY  F.  NOTE  SATURATION 
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FI6.  5  ATTENUATION  AND  TEMPORAL  DISTORTION  OF  LASER  PULSE  AS  IT 
PROPAGATES  THROUGH  SODIUM  VAPOR  CAUSING  IONIZATION 
THROUGH  RESONANCE  SATURATION.  F(t)  IS  THE  LASER  PHOTON 
FLUX  DENSITY. 


COLUMN.  LONER  LEFT  FIGURE  DISPLAYS  MEASURED  SODIUM  ATOM 
RADIAL  DISTRIBUTION  WITHIN  HEAT  SANDWICH  OVEN. 
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FIG.  7  COMPARISON  OF  EXPERIMENTAL  LASER  PULSE  AND  EMPIRICAL  FIT 
USED  IN  LIBORS  CODE. 


COLUMN  DEPTH  (cm) 


Lc**r  Pul** 


FIG.  8  SE-LIBORS  CODE  PREDICTION  OF  THE  VARIATION  IN  THE  DEGREE 
OF  IONIZATION  OF  SODIUM  VAPOR  ALONG  PATH  OF  LASERS  FOR 
THREE  INSTANTS  OF  TIME  CORRESPONDING  TO  THREE  VALUES  OF 
LASER  DETUNING. 
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FI6.  9  COMPUTED  VARIATION  OF  ELECTRON  TEMPERATURE  T.  AND 
DEGREE  OF  IONIZATION  a  ALONG  PATH  OF  LASER  PULSE  10  NS 
AFTER  START  OF  INCIDENT  IRRADIATION  FOR  THREE  SETTINGS 
OF  LASER  DETUNING. 
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FIG.  10  COMPUTED  RADIAL  VARIATION  OF  ELECTRON  TEMPERATURE  T. 

AND  ELECTRON  DENSITY  N.  FOR  10  AND  25  NS  AFTER  START 
OF  INCIDENT  IRRADIATION  FOR  A  19.7  nJ  LASER  PULSE 
(DETUNED  BY  0.14  NM)  KITH  A  GAUSSIAN  RADIAL 
DISTRIBUTION. 
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11  COMPUTED  VARIATION  OF  ABSORBED  LASER  ENERGY  FLUENCE 
WITH  INCIDENT  PEAK  PHOTON  FLUX  DENSITY  (AND  ENERGY 
FLUENCE)  FOR  A*  *  589.3  NM  AND  PEAK  SODIUM  DENSITY  OF 
2 *10 1  *  CM*1. 
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FIG.  13  SCHEMATIC  OVERVIEW  OF  NEW  Nd-YAG  LASER  DRIVEN  LIBORS 
FACILITY. 


FIG.  14  PICTURE  OF  LIBORS  FACILITY  SHOWING  HEAT  SANDWICH  OVEN 
AND  PHOTODETECTION  SYSTEM  ON  GRANITE  TABLE  IN 
FOREGROUND  AND  Nd-YAG  LASER  IN  BACKGROUND. 
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FIG.  15  SPECTRAL  SCAN  OF  EMISSION  FROM  SODIUM  PLASMA  CREATED  BY 
LASER  RESONANCE  SATURATION. 
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FIG.  16  HIGH  RESOLUTION  SPECTRAL  SCAN  OF  THE  42D-32P  MULTIPLET 
FROM  SODIUM  PLASMA  CREATED  BY  LASER  RESONANCE 
SATURATION. 


FIG.  17  LOG  PLOT  OF  EXPERIMENTAL  LINE  EMISSION  FROM  SODIUM 
PLASMA  (0).  ALSO  SHOWN  IS  BEST  LINEAR  FIT  TO  DATA  AND 
COMPUTED  LINE  EMISSION  FOR  A  ONE  SLAB  OPTICALLY  THICK 
MODEL  (A). 
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FIG.  18  EXPERIMENTALLY  MEASURED  RADIAL  DISTRIBUTION  OF  SODIUM 

EMISSION  ON  4  2D-3  2P  (_  _)  AND  62D-32P  ( _ ) 

MULTIPLETS.  [  1 


FIG.  19  DECONVOLUTED  RADIAL  ELECTRON  TEMPERTURE  DISTRIBUTION 
MEASURED  AT  60  NS  FROM  ONSET  OF  LASER  IRRADI ANCE  FROM 
CENTRE  OF  SODIUM  HEAT  SANDWICH  OVEN.  ALSO  SHOWN  ARE 
THREE  THEORETICAL  RADIAL  TEMPERATURE  DISTRIBUTIONS 
LIBORS  CODE  FOR  THREE  DIFFERENT  LASER  RADIAL 

PROFILES. 
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FIG.  20  TWO  SPECTRAL  PROFILES  OF  SODIUM  42D-32P  MULTIPLET, 
EXPERIMENTALLY  MEASURED  CLOSE  TO  THE  LASER  BEAM  AXIS 
(y  •  0.125  MM)  ANO  IN  WING  OF  PLASMA  (y  -  2.25  MM). 


FIG.  21  VARIATION  OF  EXPERIMENTALLY  MEASURED  PEAK  TRANSMITTED 
LASER  PULSE  AS  FUNCTION  OF  DETUNING  MAXIMUM  SODIUM 
DENSITY  IN  SANDWICHED  OVEN  WAS  4*10“  CM'3. 
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FI6.  22  COMPUTED  VARIATION  OF  TRANSMITTED  PEAK  LASER  PHOTON 
FLUX  DENSITY  AGAINST  LASER  SPECTRAL  DETUNING  FOR  THREE 
VALUES  OF  INCIDENT  LASER  ENERGY.  PEAK  SODIUM  DENSITY 
4*101€  CM" *. 
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590.1  NM. 
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FIG.  24  COMPUTED  RATIO  OF  ELECTRON  THERMAL 
COEFFICIENT  A  OVER  SPITZER-HARM  VALUE 
OF  ELECTRON  TEMPERATURE  Te  FOR  THREF 
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FIG.  26  NEW  SELF-PUMPING  RESONANCE  INVERSION  SPRI-CONCEPT  FOR 
SHORT  WAVELENGTH  LASER. 
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FIG.  30  COMPUTED  EMISSION  ON  THE  2.82  pM  AND  177  NM  TRANSITIONS 
OF  SrII  RESULTING  FROM  IRRADIATION  MITH  SELF-NORMALIZED 
LASER  PULSE  SHOWN.  LASER  PEAK  IRRADIANCE  WAS  1"10*W 


FIG.  31  NEW  THREE  PHOTON  SATURATION  FACILITY  SHOWING  STRONTIUM 
HEAT  CROSS  IN  BACKGROUND  AND  Nd-YAG  LASER  PUMPED  DYE 
LASER  IN  FOREGROUND. 


Physical  Mechanism  for  LORMA-Effect 


In  an  attempt  to  explain  in  a  simple  manner  the  "laser  off  resonance 
maximum  absorption"  (or  LORMA)-effect  we  shall  assume  that  a  near  dynamic 
balance  between  laser  absorption  and  superelastic  electron  quenching 
collisions  determines  the  resonance  and  ground  level  populations  of  the 
sodium  atoms.  We  further  assume  that  the  free  electron  density  Ng  is  high 
enough  that  N  >>  A^  and  we  can  neglect  spontaneous  decay  of  the 
resonance  level  populations.  Consequently,  we  can  write 


(gNj  -  n2)r21  »  Ne (N2K21  -  n1k12) 


(Al) 


An  order  of  magnitude  calculation  on  the  terms  of  the  one  dimensional 
radiative  transfer  equation  (27)  reveals  that  to  good  approximation 


dr  *  C,  -  9N,)R?,E 


i'"2r2i 


(A2) 


0 

where  I  is  the  laser  irradiance  and  the  resonance  state  energy. 
Combining  the  above  equations  leads  us  to  state  that  the  energy  extracted 
per  unit  volume  from  the  laser  pulse 


AE*  *  /  Ne(N2K21  -  N1K12)E21dt 
o 


(A3) 


l 

where  -c  is  the  duration  of  the  laser  pulse. 

In  general  equation  (Al)  can  be  written  in  the  form 


-E  /kT 

N2  g  (Sj  +  e  e) 


N 


1 


T-n- 


(A4) 


wh 


ere  Sj  =  R 21  ^NeK21  and  1S  termed  "saturation  parameter".  If  the 


]  (All) 


ti  */  \  2  rl  ^  fx  i  )  f2a-\  ^  -x7x  ,  1  -2x*x 

f (a,  x  /x)  =  azx[?  +  ( — —  )e  +  ?  e 


When  the  laser  is  only  slightly  detuned  from  resonance  x  <<  x  and  o=l. 
In  which  case 


f(a,  x*/x)  =  |  xa2  (A12) 

and  clearly  the  laser  absorption  increases  with  increased  detuning  because 

x  increases  with  increased  detuning  while  a  does  not  change  appreciably. 

0 

For  large  detuning  of  the  laser,  x  >>  x  and  a  <<  1.  Under  these 
conditions 


f(a,  x*/x)  *  (1  -  (t)ax*  (A13) 

and  it  is  apparent  that  laser  absorption  now  decreases  with  increased 
detuning  because  a  becomes  smaller  with  increased  detuning.  Also  for  large 
detuning  the  assumption  of  resonance  saturation  is  expected  to  break  down. 
It  is  thus  clear  that  there  will  exist  some  small  value  of  detuning  for 
which  the  laser  attenuation  is  a  maximum  if  electron  superelastic  collision 
quenching  plays  a  key  role  in  the  extraction  of  energy  from  the  laser 
pul se. 
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